AD  609177 


INVESTIGATIONS  ON  THE  DIRECT  CONVERSION  OF  NUCLEAR 
FISSION  ENERGY  TO  ELECTRICAL  ENERGY  IN  A  PLASMA  DIODE 


ANNUAL  REPORT 
for 

Nonr-3I09(00) 


COPY  OF  TL-g  I  ^ 


HARD  COPY 
MICROFICHE 


$.  0 

$.  /,  oo 


/I^P 


OCTOBER  31,  1964 


rvsBDC 

IjpEEUSDn 
S  «EC«3«4 

nsM  D 


RESEARCH  LABORATORIES 
GENERAL  MOTORS  CORPORATION 
WARREN,  MICHIGAN 


\ 


CLEARINGHOUSE  FOR  FEDERAL  SCIENTIFIC  AND  TECHNICAL  INFORMATION  CFSTI 

DOCUMENT  MANAGEMENT  BRANCH  410.11 


LIMITATIONS  IN  REPRODUCTION  QUALITY 


ACCESSIONS 


LEGIBILITY  OF  THIS  DOCUMENT  IS  IN  PART  UNSATISFACTORY. 
REPRODUCTION  HAS  BEEN  MADE  FROM  THE  BEST  AVAILABLE  COPY. 


Q  2.  ORIGINAL  DOCUMENT  CONTAINS  COLOR  OTHER  THAN  BLACK  AND  WHITE 
AND  IS  AVAILABLE  IN  LIMITED  SUPPLY.  AFTER  PRESENT  STOCK  IS  EX¬ 
HAUSTED,  IT  WILL  BE  AVAILABLE  IN  BLACK-AND-WHITE  ONLY. 


n  3.  THE  REPRODUCIBLE  QUALITY  OF  THIS  DOCUMENT  IS  NOT  ADEQUATE 
FOR  PUBLIC  SALE.  AVAILABLE  TO  CUSTOMERS  OF  THE  DEFENSE 
DOCUMENTATION  CENTER  ONLY. 


n  4.  DOCUMENT  AVAILABLE  FROM  CLEARINGHOUSE  ON  LOAN  ONLY 
(TECHNICAL  TRANSLATIONS). 


PROCESSOR: 


TSL-l07-i2/64 


OFFICE  OF  NAVAL  RESEARCH 
Contract  Nonr-3109(00) 


Investigations  on  the  Direct  Conversion  of  Nuclear 
Fission  Energy  to  Electrical  Energy  in  a  Plasma  Diode 

Report  No  .  5 


Authors 

C .  B .  Lef f ert 

D.  B.  Rees 

F.  E.  Jamerson 


Report  for  Period  November  1^  I963  to  October  3I,  1964 


Research  Laboratories ^  General  Motors  Corporation 

Warren^  Michigan 


Reproduction  in  whole  or  in  part  is  permitted  for  any 
purpose  of  the  United  States  Government 


This  report  has  been  prepared  under  Contract 
No.  Nonr-3109(00)  for  the  Office  of  Naval  Research. 

It  was  technically  supervised  by  Dr.  J.  J.  Connelly,  Jr., 
Commander  U .  S  .  Navy . 


ABSTRACT 


The  ionization  of  pure,  mixed  and  cesium-seeded  noble  gas  systems  by 
fission  fragments  has  been  investigated  using  ceramic-metal  tubes  operated 
in  a  high  neutron  flux  region  of  a  nuclear  reactor.  An  analysis  of  the 
production  and  transport  of  ions  and  electrons  in  the  tube  has  shown  that 
when  the  charge  loss  in  the  plasma  is  predominantly  by  volume  recombina¬ 
tion  of  molecular  ions  and  electrons,  the  tube  current  is  proportional  to 
(voltage)  .  Experimental  current -voltage  data  in  neon,  argon,  xenon 
and  a  neon-argon  mixture  with  Ar/Ne  of  10“ ^  for  gas  pressures  of  240  and 
400  torr  and  a  neutron  flux  10  ^  cm  ^  sec  ^  exhibited  the  dependence  of 
current  on  (voltage)  .  Calculation  of  the  magnitude  of  the  current  com¬ 
pletely  from  theory,  using  a  detailed  model  of  ion  generation  rate,  showed 
an  agreement  with  the  experimental  data  within  +10 %.  For  a  gas  pressure 
of  240  torr,  where  the  ion  generation  rate  was  of  order  5  x  10^^  cm’^  sec”^ 
the  average  ion  density  computed  for  the  pure  gases  was  2.^  x  10^^  cm”^ 
but  for  the  Ne-Ar  mixture  (where  the  loss  rate  of  argon  ions  is  reduced) 
the  ion  density  was  about  four  times  higher,  i  .e .  --^1.0  x  10^^  cm'^.  Some 
preliminary  studies  of  reaction  kinetics  in  Ne-Ar  mixtures,  designed  to 
maximize  ion  number  density,  have  indicated  that  with  the  kno^vn  ion  genera¬ 
tion  rates^at  a  neutron  flux  of  IC^^  cm  sec  densities  as  high  as 
^  4  X  10  cm  ^  may  be  obtained  for  a  mixture  with  Ar/Ne  of  10“^  at  values 
of  p  100  torr .  Lack  of  detailed  knowledge  of  some  of  the  relevant  con¬ 
version  processes  in  an  argon-cesium  mixture  inhibits  at  present  a  similar 
study  of  reaction  kinetics  in  this  system.  However,  preliminary  current- 
voltage  data  for  Ar/Cs  of  10  ^  at  an  argon  pressure  of  240  torr,  appears  to 
follow  a  current -( voltage  )^'^^  relationship  and  indicates  that  the  ion  gen¬ 
eration  rate  is  significantly  higher  than  that  for  Ne-Ar  or  the  pure  gases 
under  the  same  conditions  . 
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INVESITGATIONS  ON  THE  DIRECT  CONVERSION  OF 
NUCLEAR  FISSION  ENERGY  TO  ELECTRICAL  ENERGY  IN  A  PLASMA  DIODE 


INTRODUCTION 


A  noble  gas  plasma  generated  by  fission  fragment  ionization,  is  being 
investigated  for  application  in  thermionic  converters  .  The  program  has  con- 
centrated  on  an  analysis  of  the  fundamental  processes  taking  place  in  the 
plasma  generated  by  this  method  .  This  is  necessary  in  order  to  arrive  at 
the  gas  parameters  that  lead  to  optimum  plasma  conditions  for  maximizing 
electron  transport  in  the  diode  . 

In  the  previous  reporting  period/^^  data  ohuained  from  ionization 
tubes  were  presented  for  gas  fillings  of  neon,  argon,  xenon  and  neon-argon, 
along  with  an  analysis  which  described  the  I-V  characteristics  obtained  . 

In  that  report  the  fit  between  theory  and  experiment  was  made  by  adjust¬ 
ing  one  parameter,  viz.,  the  average  ion  generation  rate. 

The  present  report  describes  an  anal;y'tical  method  for  computing  the 
ion  generation  rate  directly  from  the  fission  fragment  energy  loss .  Com¬ 
bined  with  the  ion  transport  theory  developed  earlier,  this  allows  the 
calculation  of  the  I-V  characteristic  to  be  made  directly  without  any 
adjustable  parameters .  Values  of  the  ion  generation  rate  computed  for 
pure  neon  and  for  neon-argon  (Ar/Ne  of  lO'^)  are  used  in  a  detailed  anal¬ 
ysis  of  the  reaction  kinetics  of  the  neon-argon  system.  This  analysis  leads 
to  an  evaluation  of  the  electron  number  density  in  the  neon-argon  mixture 
for  varying  conditions  of  Ar/Ne  and  total  pressure .  Preliminary  data  on  ion 
generation  rate  are  also  reported  for  an  argon-cesium  mixture  with  Cs/Ar  of 
-o  3x10  .  Finally,and  incidental  to  the  main  program,  a  value  is  given  for 

the  secondary  electron  yield  per  fission  fragment  emitted  from  the  uranium 
surface,  since  that  quantity  was  readily  determined  from  the  present  work. 

The  subsequent  sections  discuss  the  program  Objectives,  Conclusions 
and  Future  Plans  . 
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OBJECTIVES 


The  objectives  initially  set  for  the  current  reporting  period  were  as 
follows : 

1.  Experiments  were  to  be  continued  on  the  inpile  measurement  of  ion 
generation  rate  in  noble  gas  mixtures . 

2 ,  Theoretical  studies  were  also  to  be  continued  on  the  ion  genera¬ 
tion  rate  and  the  ion  tube  I-V  characteristic . 

3*  Electron  transport  measurements  in  noble  gas  plasmas  were  to  be  ini¬ 
tiated  .  A  tube  incorporating  a  nuclear  heater  thermionic  emitter 
was  to  be  developed  for  this  study. 

During  this  reporting  period  considerable  time  was  spent  on  theoretical 
analyses  .  A  major  endeavor  was  the  development  of  the  ion  generation  rate 
and  I-V  characteristic  analysis.  Additionally,  a  reaction  kinetics  analysis 
was  begun  in  order  to  evaluate  electron  number  density  from  the  ion  genera¬ 
tion  rate  data  for  noble  gas  mixtures.  Consequently  the  progress  on  the 
electron  transport  experiment  was  limited  to  the  evaluation  of  a  particular 
diode  design  and  some  heat  transfer  mockup  studies  .  This  thermionic  diode 
will  be  discussed  in  a  future  report. 

CONCLUSIONS 


1.  The  current -voltage  characteristic  for  the  ceramic-metal  ionization 
tube  can  be  calculated  directly  without  the  need  for  adjustable  param¬ 
eters  .  Tlie  agreement  between  the  theoretical  and  experimental  values  of 
current  for  the  pure  gases  is  +  10 7o. 

2.  The  ion  generation  rate  can  be  computed  from  the  fission  fragment  energy 
loss  using  a  linear  energy  loss  model.  The  use  of  a  quadratic  energy 
loss  model  influences  the  resultant  I-V  characteristic  only  slightly 
when  compared  to  that  using  the  linear  model.  Computer  time  is  reduced 
by  using  the  linear  model, 

3  •  The  ion  generation  rate  can  be  computed  accurately  provided  the  value 
of  w  (ev/ion  pair)  is  Imown  for  the  gas.  For  example,  the  difference  in 
the  I-V  characteristic  (for  V  25  volts)  between  neon  and  a  mixture  of 
neon-argon  (Ar/Ne  of  10  ^)  was  accounted  for  correctly  by  the  difference 
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4. 


in  w  for  these  two  gas  compositions . 

Analysis  of  the  I-V  data  for  neon  and  neon-argon  leads  to  values  for  the 
ion  generation  rates  of  atomic  neon  ions  and  metastable  neon  atoms . 

5 .  The  electron  number  density  for  the  neon-argon  mixture  can  be  computed 
from  the  reaction  kinetic  equations  which  take  account  of  the  various 
processes  of  production  and  loss  of  ions.  The  number  density  is  sensi¬ 
tive  to  gas  pressure  and  the  ratio  of  gas  spejies.  Densities  greater 
12  -  'R  ° 

than  10  cm  were  computed  for  typical  conditions  . 

6o  Using  a  modified  ceramic-metal  ion  tube  the  ion  generation  rate  in  an 
argon-cesium  mixture  (Cs/Ar  of  3  x  10  was  shown  to  be  higher  than 
for  pure  Ar  and  for  Ne-Ar, 

7*  Prom  a  reevaluation  of  data  taken  with  a  vacuum  diode  operated  inplle^ 
the  secondary  electron  yield  from  the  uranium-nickel  (5.7  wt.^o  Ni)  alloy 
used  in  the  ion  generation  rate  studies  is  207  +  10  electrons  per  fission 
fragment  emitted  from  the  foil. 

8 .  The  ion  generation  rates  and  plasma  densities  at  the  midpoint  of  the  ion 
tube  (3  mm  gap)  for  a  gas  pressure  of  240  torr  and  a  neutron  flux  of 
1.0  X  10  ^  cm  ^  sec  ^  are  given  below. 


Gas 


Neon 


Arg'ju 


Ion  Generation  Rate  Recombination 

(cm  sec  ) _ _ Loss  Process 

1 .8  X  10^^ 


3.6  X  10^^ 


y  +  24  — ^  y^  +  y 
^  +  e“ - >2y 


Ion  Density 
(cm"^) 

2.9  X  10^^ 


2.3  X  10^^ 


Xenon 


8.1  X  10^^ 


2.4  X  10^^ 


Neon-Argon  2.6  x  10 
(0.1%  Ar) 


INe  as  above 

Ar”^  +  Ar  +  Ne — +  Ne  1 .0  x  10^^ 
Ar^^  +  e  — >  2  Ar 


The  results  for  S  are  higher  in  this  case  than  before^^^  and  are  attri¬ 
buted  to  a  more  precise  evaluation  of  the  tube  geometry  in  calculating  the 
ion  generation  rate .  The  value  of  number  density  for  neon-argon  was  computed 
from  the  reaction  kinetics  analysis  . 
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FUTUEE  PLANS 


1.  Rsaction  kinetics  studies  will  be  extended  to  include  an  analysis  of  the 
argon-cesium  plasma  ion  number  density. 

A  microwave  measurement  of  ion  number  density  in  fission  fragment  ionized 
mixed  gas  systems  (Ne-Ar  and  Ar-Cs)  will  be  made  inpile.  This  experi¬ 
ment  wii  1  be  used  to  verify  the  reaction  kinetics  analysis  . 

3-  A  diode  using  a  nuclear  thennionic  emitter  with  a  mixed  gas  filling  (near 
optimum  conditions)  will  be  operated  inpile  in  order  to  determine  the 
electron  transport  through  a  noble  gas  plasma. 

4 .  Analysis  of  electron  transport  through  noble  gas  plasmas  related  to  the 
inpile  electron  transport  measurement  will  be  initiated . 

REFERENCES 
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Electrical  Energy  in  a  Plasma  Diode,  I963  Annual  Report  No.  4;  C.  B. 
Leffert,  F.  E.  Jamerson,  and  D.  B.  Rees. 

ACKNOWLEDGMENTS 

The  authors  are  indebted  to  Professor  D.  J.  Rose  of  Massachusetts 
Institute  of  Technology  for  his  stimulating  discussions  and  his  formula¬ 
tion  of  the  reaction  kinetics  problem  for  the  mixed  gases. 

The  assistance  of  the  Universitv  of  Michigan  reactor  staff  during  the 
course  of  these  experiments  is  gratefully  acknowledged. 

Technical  assistance  in  the  design  and  fabrication  of  the  tubes  and 
inpile  test  system  has  been  given  by  Messrs.  R.  Aikin,  A.  Dolenga,  R.  Dus- 
man,  R.  Knoll,  D.  Lee  and  J.  Palazzolo. 


k 


BLANK  PAGE 


SECTION  A 


ION  GENERATION  RATE  THEORY 
AND  ANALYSIS  OF  EXPERIMENTAL  CURRENT- VOLTAGE  DATA 

ABSTRACT 

The  ion  generation  rate  in  a  noble  gas  is  computed  from  the  energy- 
loss  of  the  fission  fragments  as  they  penetrate  the  gas  adjacent  to  a 
fissioning  uranium  plate.  A  description  of  the  digital  computer  pro¬ 
gram  is  given.  The  current -voltage  characteristic  for  an  ionization 
tube  is  then  computed  using  the  ion  generation  rate .  A  brief  descrip¬ 
tion  of  this  computer  program  is  also  given.  Finally,  the  inpile 
experimental  current -voltage  characteristics  obtained  on  a  metal - 
ceramic  ionization  tube  are  compared  with  the  results  of  the  compu¬ 
tations  .  The  I-V  curves  computed  with  no  adjustable  parameters  agreed 
well  (+10%)  with  the  experimental  curves  and  confirmed  both  the  ion 
generation  rate  and  ion  transport  models . 
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I .  INTRODUCTION 


The  fission  fragment  ionization  of  noble  gases  has  been  under  study  for 
some  time  at  the  General  Motors  Research  Laboratories  as  a  means  to  create  a 
space  charge  neutralizing  plasma  in  a  thermionic  energy  converter .  Current - 
voltage  characteristics  have  been  obtained  on  ionization  tubes  run  in  the 
University  of  Michigan  nuclear  reactor  and  from  those  data  it  was  possible  to 
deduce  an  average  value  of  the  ion  generation  rate  in  the  fission  fragment 
generated  plasma .  ^  It  was  realized  early  in  these  studies  that  it  would  be 
very  helpful  if  the  ion  generation  rate  at  a  point  in  the  gas  could  also  be 
obtained  completely  from  theory  starting  v/ith  the  energy  loss  of  the  fission 
fragments  and  taking  into  account  the  actual  geometry  of  the  tube .  As  a  first 
step  towards  this  goal  an  analytic  solution  for  the  ion  generation  rate  was 
ODtained  for  the  ideal  case  of  infinite  parallel  plane  geometry  and  the  solu¬ 
tion  was  discussed  previously 

This  report  describes  the  computation  of  the  ion  generation  rate  for  the 
two  dimensional  (R,Z)  geometry  of  the  ionization  tube  and  the  digital  computer 
code  employed  in  the  solution.  In  order  to  compare  theory  with  experiment  a 
second  program  was  written  to  compute  the  current -voltage  characteristic  for 
the  ionization  tube.  This  second  code  employs  a  generalized  treatment  of  the 
ion  transport  model  ^  and  takes  as  input  the  computed  ion  generation  rate 
from  the  first  code  for  a  matrix  of  points  in  the  interelectrode  gap  and  then 
computes  the  sheath  thickness  and  current  for  specified  voltage  points  .  The 
program^details  for  the  first  versions  of  these  codes  are  presented  in  another 
report  and  only  the  modifications  for  the  final  version  are  described  here . 

Finally  the  fit  of  the  theory  to  experiment  is  discussed  for  the  pure 
noble  gases  and  for  the  mixture  of  neon  and  argon . 

II  •  ION  GENERATION  RATE  FROM  FISSION  FRAGMENT  ENERGY  LOGS 

The  energy  loss  of  fission  fragments  as  they  penetrate  matter  has  been 
treated  theoretically  by  Bethe  and  Ashkln.^^^  The  fission  fragment  loses 
energy  rapidly  by  electronic  collisions  at  the  beginning  of  its  track  when  its 
velocity  and  charge  are  greatest .  As  the  fragment  slows  dora  it  loses  part 
of  its  charge  through  recombination  v/ith  electrons  and  the  energy  deposition 
rate  decreases.  Nuclear  collisions  become  significant  near  the  end  of  the  track; 
this  increases  the  energy  deposition  rate  but  not  necessarily  the  ionization  rate. 
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Ionization  of  Gases 


The  ionization  of  gases  by  fission  fragments  has  been  measured ^ ^ and 

these  results  have  been  analyzed  in  terms  of  an  "ionization  defect"  A 

i.e.  E^=wI-rA  where  is  the  initial  fragment  energy^  I  is  the  number  of  ion 

pairs  formed  by  fragments  stopped  in  the  gas  and  w  is  the  energy  loss  per  ion 

pair  for  these  high  velocity  fragments.  The  ionization  defect, A  ,  accounts 

for  the  non-ionizing  energy  lost  in  nuclear  recoils  and  w  has  been  found  to  be 

very  nearly  equal  to  the  value  w  =E  /l  obtained  for  alpha  particles  of  initial 

energy  E^  in  the  same  medium.  A  similar  quantity  for  the  fission  fragments 

is  w=E^/l  and  setting  w'^w^  we  have  w/w^=l+ A /iw^.  The  fission  fragments  are 

divided  into  two  median  energy  groups  (light,  E  =98  Mev  and  heavy,  E  =67  Mev) 

01  o2 

and  the  average  energy  loss  per  ion  pair  w  relative  to  v  is  w  =1.06  w  for  the 

(h\^  I  a 

light  fragment  and  W2=l  .11  w^  for  the  heavy  fragment.''  ^ 

The  ionization  rate  along  the  fission  fragment  track,  dl/dx,  can  be  writ¬ 
ten  in  terms  of  the  stopping  power  of  the  gas  ( -dE/dj^  as  (dl/dx)  .=(-dE/dx)  .x 
(dl/dE)j  where  the  subscript  j  equals  1  and  2  for  the  light  and  heavy  energy 
groups  of  fission  fragments,  respectively.  In  this  study  the  quantity  w  is 
made  independent  of  E  over  the  entire  energy  range  of  the  fragment,  (dl/dx) .= 
(-dE/dx)j/wj,  and  so  the  ionization  rate  is  expressed  in  terms  of  the  stopping 
power  and  the  eneigy  loss  per  ion  pair  for  cc  particles  in  the  same  gas. 


Energy  Loss  Model  for  Fission  Fragments 

The  energy  of  the  fragment  (E  )  as  a  function  of  distance  x  from  birth  has 

(7) 

been  approximated  by  either  a  linear  dependence  on  distance,  E  =E  (1-x/R  ), 

p  j  oj" 

or  a  quadratic  dependence  on  distance  E^=E^^ (l-x/R.  .)  here  R.  ,  is  the  total 

range  of  the  fragment  j  in  the  medium  i  .  In  terms  of  range  units ,  u  =x  /R 

ij  i'  ij^ 

the  energy  deposition  per  range  unit  for  the  linear  model  is 


dE^/da4,.  =  (1) 

and  is  independent  of  the  previous  history  of  the  particle,  while  the  energy 
deposition  per  range  unit  for  the  quadratic  model  is 

dFi/du,.  =  -2eo^(i -5  (2) 

and  depends  upon  the  distance  traveled  in  the  f  lel  (iti)  as  well  as  in  the  gas 
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(i=2)..  For  the  conditions  of  interest  here,  that  is  for  fuel  thicknesses  greater 
than  one  fission  fragment  range  in  the  uranium  and  for  gas  film  thicknesses  con¬ 
siderably  less  than  one  fission  fragment  range  in  the  gas,  it  is  not  expected 
that  the  final  ion  generation  rate  \^d.ll  be  very  sensitive  to  the  particular 
energy  loss  model  employed  since  at  any  one  point  in  the  gap  the  gas  "sees" 
fragments  from  nearly  all  portions  of  the  fission  fragment  energy  spectrum. 

Range  of  Fission  Fragments 

The  fission  fragm.ent  range  is  inversely  proportional  to  the  density 
p(gm/cm  )  of  the  stopping  medium. For  a  wide  range  of  materials,  the  stop¬ 
ping  power  was  found  to  be  proportional  to  the  square  root  of  the  average  atomic 
mass  A  of  the  stopping  medium^^^  to  give  the  following  relationship^^^  for  the 
unknown  range  R  in  terms  of  the  known  range  R^  in  air:  R=R  (p  /p)(A/A  )^'^^cm. 

For  a  gas  at  temperature  T(°K)  and  pressure  p(torr),  this  reduces  to 
R=20.1  Rg_(  V^/M)(T/p)cm,  where  M  is  the  molecular  weight  of  the  gas. 

Derivation  of  Ion  Generation  Rate  Equation 


Fig.  1.  Geometry  for  derivation  of  ion  generation  equation. 
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Tlie  ion  generation  rate  Sir^,z^)  at  point  P  in  the  gas  (see  Pig.  l)  from 

both  fission  fragments  is  obtained  from  the  energy  deposition  rate  £.  {eV ,cm~^sec~^ 
at  P  and  w.,  the  energy  to  create  an  ion  pair  in  the  gas 

a. 

SCk^Zz.)  -  2  .  (3) 

Consider  the  differential  energy  deposition  rate  at  P  resulting  from 

those  fragments  of  group  j  that  originate  in  the  volume  element  dv^^  with  coor¬ 
dinates  r^,\\J  in  the  uranium.  The  fission  rate  in  dv^  is  (Q^Mv^)  where  Q  is 
the  fission  cross  section  (cm  and  4)  is  the  thermal  neutron  flux  (cm"^sec"^)^ 
and  the  total  kinetic  energy  released  in  dv^  to  fragments  of  group  j  is 
(Qf <t)EojdVi)  where  is  the  initial  fragment  energy.  Assuming  isotropic  emis¬ 
sion  of  the  fragments  within  the  uranium  and  straight  line  flight  over  the 
entire  range,  the  fraction  of  fragments  that  would  penetrate  dv^  at  P  with  no 
absorption  is  dA2COs^O/l4.jr(z^+z^)^ .  The  energy  deposition  per  cm  of  path  at  P 
for  each  fragment  from  dv^  is  given  by  Eqs .  (l)  or  (2),  depending  upon  the 

energy  loss  model,  and  the  fraction  of  the  initial  kinetic  energy  deposited  in 

dVp  is  (dE  /du  )(dz  /R  E  .cosO) . 

^  J  2  2j  oj 

Combining  these  factors,  the  differenti.al  ion  generation  rate  at  P  due  to 
the  fission  fragments  of  group  j  from  the  volume  element  dv^  in  the  fuel  is 


dS-(z,,z._e) 


F-  ,-u.')de 


where  for  the  linear  energy  loss  model 


fz*  (e,  J  =  tao  Q 

and  for  the  quadratic  energy  loss  model 


(5) 


=  2  //  -  3^'  ±Jd‘i.  ]t^n&  . 

V  C<K>©  J 

For  an  infinite  slab  of  uranium  of  thickness  tS:  R^^  the  integrations 
over  the  coordinates  z,  and  0  are  readily  performed  analytically  to  give 
solutions  to  Eq .  (4): 


6 


Linear  Energy  Loss^^^ 

~L  ■  (,  ) 


[t  (|  -<  (8) 

A  comparison  of  the  two  energy  loss  models  for  this  ideal  geometry  is  given 
in  Fig,  2.  Since  the  constant  is  the  same  in  Eqs .  (7)  and  (8)  only  the  factors  in 
brackets  (-f  in  Fig.2)ai‘e  plotted  as  well  as  their  ratio  for  a  gas  pressure  of 
2^1-0  torr .  Various  dimensions  in  the  tube  are  show  on  the  abscissa  for  various 
gases.  If  an  average  dimension  is  the  interelectrode  spacing,  d,  there  is  less 
than  20°/c  difference  in  ion  generation  rate  predicted  between  the  two  models. 

For  a  point  P  in  the  gas  off-axis  and  above  a  finite  uranium  disc 
the  integrations  over  the  coordinates  z,Vi/  and  Q-  become  quite  complex.  The  com¬ 
plexity  arises  from  the  geometry  and  finite  range  of  the  fission  fragments  and 
appears  in  the  limits  of  integration.  The  solution  to  Eq.(4)  takes  the  foirni 


^  J 

Quadratic  Energy  Toss 

=  E  F-i') 

\  wj  7 


It  Jr 


where  some  of  the  limits  of  integration  change  depending  on  whether  ^Si’2-S.Pq 
or  Pp  — ^2“ ^2’  ^  defined  by  Eqs.(5)  or  (6)  depending  upon  the 

energy  loss  model.  The  subscript  k  designates  a  range  of  integration  over  the 
variables  u^^  and  Vi/ and  this  is  discussed  in  the  next  section.  The  expression 

to  be  used  for  the  above  limits  of  integration  for  each  value  of  k  are  listed 
in  Ta,ble  1 . 
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mBLE  1 


Limit 

k  =  1 

Value  Eq . 

k  =  2 

Value  Eq . 

k  =  3 

Value  Eq . 

k  =  4 

Value  Eq . 

Case  for  0:s.r2S.P^ 

min 

.. 

Ijk 

0 

“5  (12) 

U“J  (12) 

(14) 

max 

""ijk 

u”J  (12) 

u”2  (1^) 

u”^  (14) 

"!j  (15) 

min 

'^■v  = 

Jlc 

0 

0 

max 

'V,  (13) 

J 

0 

max 

Jk 

Jt 

(13) 

J 

IT 

It 

^min 

jk 

0 

0 

0 

0 

^max 

jk 

^max 

Pi 

gmax 

Pi 

(10) 

J 

Case  for 

min 

Ijk 

0 

(12) 

uf.  (12) 

(18) 

max 

""ijk 

(18) 

(IS) 

“!j  (14) 

min 

H/jk  = 

(17) 

(17) 

vy7^  (13) 

vf" 

H/r  = 

n 

\y7^  (13) 

n 

It 

^min 

ok 

^min 

P^  (l6) 

(16) 

Pi 

(16) 

Pi 

PIHH 

^max 

jk 

^max 

(11) 

. . .  1....  <„ 

^max 

Pi 

(10) 

J 

^max 

J 
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Limits  of  Intep, ration 

The  detailed  derivation  of  the  limits  of  integration  will  not  be  given 
here  but  only  an  outline  of  the  derivation  together  with  the  final  expressions  . 
The  limits  of  integration  are  fixed  by  the  finite  range  of  the  fission  frag¬ 
ments  .  This  range  together  with  the  geometry  iafines  the  volume  of  integra¬ 
tion  in  the  uranium  that  can  contribute  to  the  ionization  at  P .  (Fig.  l)  For 
a  given  value  of  and  Vj/  the  integration  over  Q-  is  performed  first  followed 
by  the  integration  ovei  and  finally  by  the  integration  over  z^ .  Because  of 
symmetry  the  upper  limit  on  Vis  taken  as  n  and  the  result  is  multiplied  by  2. 

It  is  helpful  to  consider  first  the  upper  limit  on  O  which  prevails  when 
O  is  not  limited  by  the  finite  size  of  the  uranium  disc.  For  any  particular 
point  in  the  gas  a  distance  z^  above  the  uranium,  fission  fragments  passing 

through^gis  point  could  originate  from  a  depth  z^  in  the  uranium  out  to  an 
angle  O,  given  by 


(10) 


Equation  (lO)  describes  a  circle  on  a  plane  at  depth  z^  in  this  uranium  of 
large  extent  and  many  of  the  limits  of  integration  are  defined  by  the  inter¬ 
action  of  this  circle  and  the  edge  of  the  actual  uranium  disc . 

Case  for  r^^ 


Fig.  3.  Geometry  for  limits  of  integration; 


case  for  0  s  r^  < 
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Referxing  to  Pig.  3,  consi'^or  the  point  F  in  the  gas  a  distance  from  the 
axis  of  the  tube  but  still  over  the  uranium  disc,  i  .e  .  r^  <  and  at  a  height 
Zg  above  the  uranium  where  z  <Rp..  For  the  first  region  of  the  uranium 

°  ^  (k=l  in  Eq.(9)),  O.  >ej^  in  Fig.  3  and  the  upper 

limit  on  O  is  fixed  by  the  edge  of  the  uranim  disc  and  is  given  by  o’^®'^(see 
Fig.  3)  where  ^1 


G 


PI 


The  range  of  integration  on  VJ/ls  0  2  \)y  -  „  and  the  limit  is  fixed  by  the 
condition  that  O*?  Oj,  which  yields  ^ 

u 


For  the  second  region  of  the  uranium  (u‘'^  <  u^^  ^  u^p  the  circle  described 
Ey  Eq.  (lO)  intersects  the  uranium  disc  and  two  regions  of  integration  over  VjV 
must  be  considered.  For  the  first  Vjf  region  (k=2  in  Eq.(9)),  i  .e .  0  ^  \\J^, 

the  range  of  integration  for  O  is  0  <0  <  while  for  the  second  \|/  region 

(k=3  in  Eq^.(9))^  i  .e .  .  <  V(y<jt^  the  range  of  integration  for  O  is 


max 


0  <  ^  <  ^max 

j 

and  solved  for  to  give 


To  obtain  (Eq.  (n))  is  set  equal  to  (Eq.(lO)) 


^1 


J 


C^''f  1 


.'(13) 


The  upper  limit  on  u^^  for  this  range  is  found  by  setting  0™^^=  0  in  Fig .  3 
which  yields 


(14) 
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For  the  third  region  of  the  uranium  (u™^<^u^^<u™3)^  k=4  in  Eq.(9)  and  the 
circle  defined  by  Eq.(l0''  lies  within  the  uranii.im  disc.  For  this  region  the 
range  of  integration  on  0  is  and  the  range  of  integration  on\\J  is 

^  upper  limit  u^^  is  the  difference  between  1  range  and  the  dis¬ 

tance  to  P  in  range  units ^  i  .e  . 


I  --u.. 


Case  for  <  r^ 


(15) 


Fig.  4.  Geometry  for  limits  of  integration;  case  for  p^  r^  <  p 

Referring  to  Fig .  4  for  r^  =»  p^  it  is  seen  that  the  lower  limits  for  © 
andVJ/  are  greater  than  zero.  For  the  first  region  of  the  uranium  O^su 
when  ©^  =.©^  and  k=l  in  Eq.(9),  the  range  of  integration  on  ©  is 

Equation  (ll)  still  defines  ©^^^  at  the  far  side  of  the  uranium  disc  for  angle^ 
V|/  and  ©p^  is  set  by  the  near  side  of  the  uranium  disc  for  angle  \^ ,  i  .e . 


e 


Vnm 


COfO 


(16) 
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The  range  of  integration  on  V|/  is  where  is  fixed  by  the  tan¬ 

gent  to  the  edge  of  the  disc  from  the  projected  point  i  .e . 


-  I 


j- 


(17) 


The  upper  limit  is  given  by  Eq.  (l2) . 

The  second  region  of  the  uranium,  u^<u^  is  now  bounded  by  a  dif¬ 

ferent  upper  limit.  Again  the  circle  described  by  Eq .  (lO)  intersects  the  disc 
and  gives  two  regions  of  integration  over  Vjy.  For  the  first  Vj;  region  (k=2  in 
El-  (9)):,  i.e.\|;  '  <  V|y<  the  range  of  integration  for  O  is 

^TniiT  ^  Tnin. 

pi  “■  ~  pi  ^p2_  given  by  Eqs .  (16)  and  (ll)^  respectively, 

Por  the  second  V|/ region  (k=3  in  Eq.  (9)),  i  .e .  the  range  of  inte¬ 
gration  for  «  Is  where  is  given  by  Eq.  (lo)  .  The  limits 

^  j  andVj;^  are  given  by  Eqs .  (l6)  and  (13),  respectively.  Tlie  limit  u 
is  given^by  Eq.  (l2)  .  The  limit  u^^  is  obtained  by  setting  for 

and  solving  for  u.,  ,  to  give  ^  ^ 

J  Ij 


For  the  third  region  of  the  uranium  (k=4  in  Eq.  (9)),  u^^^u  the 

range  of  integration  on\|./ is  m  n  and  the  range  of  integration 'on  ^ 

min  max  J  ' 

^pP  ~  whei e  all  of  these  limits  have  been  defined  previously.  A 

summary  of  the  limits  of  integration  for  Eq,  (9)  is  given  in  Table  1. 


Computer  Programs  and  So.lutions  for  S(r2,Z2) 

The  solution  of  Eq .  (9)  for  S(r2,Z2)  requires  numerical  techniques.  Two 
programs  were  written  for  solution  on  the  IE  I  709^  computer.  Tlie  first  pro¬ 
gram  for  linear  energy  loss  model  is  described  in  detail  in  another 

report.  The  second  was  a  more  general  program  and  was  used  for  the  quad¬ 

ratic  energy  loss  model.  Both  versions  are  described  briefly  together  with 
some  representative  output . 

All  numerical  Integrations  vrere  performed  by  Simpson's  rule,  using  a 
fixed  number  of  steps  (input)  for  each  vai’iable  .  Since  the  limits  of  inte¬ 
gration  were  variable,  the  increment  size  changed  from  integration  to 
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integration.  One  subroutine  was  written  which  performed  integrations  for  both 
analytic  and  tabular  functions  having  been  given  the  name  of  a  one  argument 
Fortran  function  as  part  of  its  argument  list .  When  the  function  was  tabular^ 
two  different  Fortran  i unctions  could  be  used  .  The  first  function  was  used 
when  the  range  (and  thus,  the  increment)  for  the  integration  was  known  in 
advance  .  The  second  function  used  an  interpolation  routine  and  the  vectors 
for  the  independent  variable  X  and  dependent  variable  F(x)  were  specified  in 
the  control  array.  A  call  to  the  function  with  argument  A  caused  the  inter¬ 
polation  routine  to  interpolate  the  argument  A  into  the  sequence  X  and  return 
the  interpolated  value  F(a)  . 

Linear  Energy  Loss  Model  For  the  linear  energy  loss  model,  F.=tan  Q- 
(see  Eq .  (5))^  and  Eq.  (9)  can  be  simplified  by  some  analytical  integrations. 
Equation  (9)  cati’ be  written  as 


4  '''i  Rh  f  ■  ^  =  I;  ^  \ ,  y.^de  . 

y.r'1  J 


(20) 


The  quantities  on  the  right  hand  side  of  Eq„  (20)  have  already  been  defined. 


After  some  analytical  integration?  Eq.  (20)  rcJuced  to  the  follo\d.ng 


expres¬ 


sions  for  f  . : 

^  J 


For  case  0^r^«c 


( dZ, 


rr> 

-^'3 


y  j 


^ J  -Iyi  Q  ^  ij/)  d  ijj~^  dZj 


7,  ’  y.j 


jy-  (5> 


+  ZTr  ^1  -  f  V^'  i-  v7']( /  -  A  ^  2/”.’ 


where 


;  B  6^/,  It-,  =  CmV  0^ 

for  case  p^ 


7,^0 


^*1  ]^'3 

Ar  /i'. 


74j  fi4.; 


o 

^'3  TT 


V-^- 


df  6l 


j  ^^1  j  ^ 2/ ,  ^  ^  ij/ 

■1  O  T  i  y  /HZ'! 

7,  =  27, .  -y,  (2^. 


where 


^  CW  (9, 


(2I4.) 
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(25) 


+  CjT^tp - 


# 


ThcsG  Gqiiations  ware  solvGd  numGrlcally  with  the  Fortran  code  and  the  inputs 
output  and  flow  diagrams  are  discussed  elsewhere. ^ 

Examples  of  the  output  of  this  code  are  shown  in  Figs .  5  and  6  .  A  linear 

plot  of  the  variation  of  S  with  the  distance  z,^  from  the  uranium  plate  is 

shown  in  Fig .  5  at  various  distances  r^  from  the  center  of  the  uranium  disc . 

Typical  values  of  the  tube  parameters  were  used,  i  .e .  xenon  gas  at  240  torr 

11  -2  -1 

pressure  and  <('  =  1.29x10  cm  sec  .  Solutions  out  to  one  full  range  of  the 
fission  fragment  in  the  gas  are  shown  in  the  log-log  plot  in  Fig,  6.  In  both 
figures  the  code  solutions  are  the  solid  curves  and  the  points  were  obtained 
from  analytic  expressions  presented  earlier  in  Section  C  of  reference  (l) . 

The  analytic  solution  along  the  axis  of  the  tube  indicated  by  full  circles 
Tivas  used  to  verify  the  accuracy  of  the  code  solution.  The  analytic  solution 
for  the  edge  of  the  uranium  disc  {v^=p^^  ,9^2)  involved  an  approximation^^ ^ 
and  the  agreement  with  the  code  is  considered  satisfactory. 

The  variation  of  the  ion  generation  rate  for  the  experimental  tube  geom¬ 
etry  (described  later)  is  shown  in  Fig.  7  for  the  case  of  neon  gas  at  a 
pressure  of  240  torr  and  a  neutron  flux  of  1  .OxlO^^m'^sec"^ .  Here 
is  plotted  versus  r^  for  three  values  of  z^,  i  .e .  Zg  =  0,  d/2  and  d.  This 
plot  shows  that  near  the  electrode  surface  (z^  =  O)  but  past  the  edge  of  the 
uranium  disc  (r^  the  ion  generation  rate  drops  very  rapidly  because 

fission  fragments  emitted  at  O  =  it/2  are  absorbed.  The  ion  generation  rate 
along  the  opposite  electrode  is  reduced  because  of  absorption  of  the  fission 
fragments  in  the  gas  but  S  is  somewhat  more  uniform  since  the  fragments  can 
now  reach  to  the  outer  edges  of  the  tube.  Values  for  the  constants  used  in 
these  computations  are  listed  in  Table  2  in  Section  IV. 
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Ion  generation  rate  -  S(r..Z.)  X  10  ci.i'^sec 


Distance  from  uranium  surface.  (cm) 


Pig .  5  •  Spatial  variation  of  ion  generation  rate  for  t:/pical 
tube  parameters :  linear  plot . 


generation  rate  -  S(r«,  Z«)  cm  "sec 


Variation  of  ion  generation  rate  for  experimental  tube  geometry 


Generalized  Energy  Loss  Model  The  solution  of  Eq.  (9)  using  Eq .  (6)  for 
the  quadratic  energy  loss  model  required  a  triple  integration  using  numerical 
techniques.  A  new  program  for  the  function  S{r^,z^)  was  coded  as  a  subrou¬ 
tine  and  substituted  into  the  main  control  program  for  the  linear  energy  loss 
subroutine  for  S{r^,z^) .  The  new  subroutine  for  was  written  for  a 

generalized  energy  loss  model  by  making  F  (G,u  u  )  a  Fortran  II  function 
statement .  ^ 

A  solution  for  S^r^^z^)  using  this  code  and  the  quadratic  energy  loss 
model  from  Eq.  (6)  is  shown  in  Fig.  5  as  the  dashed  curve  for  S  along  the 
axis  of  the  tube  .  The  ion  generation  rate  is  seen  to  be  the  same  at  the 
uranium  surface  (see  Eqs  .  (7)  and  (8)  in  limit  u^j^O)  and  do\^m  only  about 
10 /o  at  the  far  electrode  (z2=d).  As  a  check  on  the  accuracy  of  the  triple 
integration  subroutine^  Eq.  (5)  was  substituted  for  the  Fortran  function 
statement  ^  curve  was  obtained  identical  to  the  solid 

curve  labeled  '^2^  .001  for  the  first  linear  energy  loss  model  code , 

Some  trouble  was  experienced  with  the  numerical  solution  of  the  integral 
in  Eq.  (9)  for  positions  in  the  gas  very  close  to  the  uranium  surface  (z^ — >0). 
This  condition  involves  the  integration  of  the  function  tan  O  over  the  angle  O 
to  an  upper  limit  approaching  jt/2  where  tan  0 — >co.  The  Simpson's  rule  formu¬ 
lation  is  a  poor  approximation  for  this  condition  and  the  problem  was  partially 
circumvented  by  introducing  a  cut-off  to  the  upper  limit  of  G.  This  problem 
did  not  arise  in  the  single  integral  solution  for  the  linear  energy  loss  model 
because  the  integration  over  0  was  done  analytically. 

It  will  be  shown  later  in  Section  III  that  the  experimental  verification 
of  these  computations  of  the  ion  generation  rate  involves  the  measurement  of 
an  ion  current  drawn  from  either  of  two  sheaths  separating  the  plasma  from 
the  two  bounding  metal  electrodes  .  When  the  sheath  is  at  the  electrode  oppo¬ 
site  the  uranium  this  problem  of  the  numerical  integration  over  tan  G  does 
not  occur  and  verification  of  the  theory  can  be  obtained  .  The  running  time 
on  the  computer  was  much  greater  for  the  triple  integral  solution  than  for 
the  single  integral  solution  for  the  linear  energy  loss  model  and  since  the 
difference  in  answers  was  less  than  10 no  further  effort  was  made  to  correct 
the  triple  integral  solution  for  z^ — ^0. 
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Ill .  ION  TRANSPORT  CURRENT-VOLTAGE  CHARACTERISTIC 


In  the  previous  section  it  was  shovm  that  the  ion  generation  rate  in  a 
gas  film  ad.iacent  to  a  fissioning  uranium  plate  could  be  computed  from  known 
parameters  such  as  the  range  of  the  fission  fragments  and  the  energy  to  create 
one  ion  pair  in  the  gas  .  It  was  of  interest  to  verify  this  theory  with  an 
experimental  measurement .  The  usual  practice  to  measure  ion  production  in 
typical  ionization  tubes  is  to  apply  a  strong  electric  field  betvreen  two  elec¬ 
trodes  and  collect  the  ions  on  one  plate  and  the  electrons  on  the  other  as 
rapidly  as  they  are  made  .  The  plasmas  of  interest  in  these  studies  are  much 
too  dense  to  bring  about  such  complete  charge  separation  so  it  was  necessary 
to  develop  a  theory  for  ion  transport  in  the  sheath  bounding  a  dense  plasma . 

Ion  Transport  Model  for  Uniform  Ion  Generation  Rate 

The  ion  transport  model  for  the  case  of  a  uniform  ion  generation  rate  (s) 
in  the  interelectrode  volume  has  been  described  previously 


Fig.  8.  Potential^  ion  density  and  field  distribution  for  the  plasma  model. 
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Briefly^  the  predominant  ion  loss  process  in  the  bulk  of  the  plasma  (see 
Fig.  8)  was  considered  to  be  dissociative  volume  recombination  so  that  for  the 
pure  gases  the  ion  density  was  given  by 

n  -  -  n"  -  (sMr) 

vz-here  CC^is  the  dissociative  recombination  coefficient.  The  applied  voltage  was 
assumed  to  appear  across  the  cathode  sheath  in  order  to  drag  the  positive  ions 
out  of  this  sheath  and  onto  the  negative  electrode.  It  was  assumed  th&L  vol¬ 
ume  recombination  could  be  neglected  in  the  sheath  as  well  as  diffusion  and 
the  basic  particle  transport  equations  in  the  region  reduced  to 


r  ~  =  - yU'  t  n~ 

and 


(27) 

(28) 


where  T  is  the  current  density  of  particles  and  p  the  mobility  coefficient. 
Using  Poisson’s  equation  and  the  continuity  equations  the  total  current  I  was 
derived  in  terms  of  the  applied  voltage  V  as 


X  =  /)  5^^* 


where  A  is  the  area  of  the  negative  electrode,  is  the  permittivity  of  vacuum 
and  e  is  the  electron  charge.  Of  particular  Importance  was  the  dependence  of 
the  current  on  the  uniform  ion  generation  rate  to  the  3/4  power  and  the  depend¬ 
ence  on  the  voltage  to  the  l/2  powex  . 


Ion  Transport  Model  for  Non-Uniform  Ion  Generation  Rate 

For  this  iOii  transport  model  the  ion  generation  rate  S(r2,Z2)  is  used. 

It  is  expected  that  the  sheath  thickness  will  be  relatively  uniform  over  the 
electrode  surface  so  that  radial  diffusion  of  the  ions  can  be  neglected.  In 
this  case  the  ion  transport  equations  can  be  solved  as  a  one  dimensional  prob¬ 
lem  and  the  distribution  of  potential,  ion  density  and  field  between  the  two 
electrodes  of  the  ion  tube  can  be  assumed  to  be  similar  to  the  sketches  of 
Fig .  8 . 
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Tne  basic  particle  transport  equations  are  Eqc  .  (27)  and  (28)  and  the  con- 
tinuity  equations  are 


dr 

d  z. 


=  dJ2^] 

dZ  /k 


S  (^1  >Zi.)  . 


(30) 


Equations  (28)  and  (30)  when  integrated  yield 


E 


SU')^ 


(31) 


which  evidently  can  satisfy  the  approximate  boundary  conditions  E  — >E  ^  0 
when  z — $>0  for  finite  n'^ .  For  b«=-  z«=d,n^:.'s(z)/Q:J^/^  .  ^ 

Since  n  »  n  in  the  sheath  Poisson's  equation  becomes 


6E 


e  n  t- 

combining  Eqs .  (31)  and  (32)  with  the  boundary  conditions  E=0  when  z=b^  w 
find  the  electric  field  distribution  in  the  sheath 


(32) 


le 


E{z)  =. 


dz'  /  S(^^")dz" 

V  'O 


/z 


'/z. 


(33) 


■oy^ 


Ihe  potential  drop  V, 
once  more^  to  give 

r  ^(>^0 


3'Cross  the  cathode  sheath  is  found  by  integrating 


v; 


b(v-i.) 


d  Z 

o  i—  o 


4-  y'  ^ 

fciz' 

L  o 


(  .) 


which  is  very  nearly  equal  to  the  potential  drop  V  across  the  ionization  tobe 
for  b(r2)  <d.  Substituting  V  for  ^  and  using  the  ion  generation  rate 

Pi’om  Eq .  (9)»  we  solve  the  integral  Eq,  (34)  to  obtain  b(r2). 

The  total  ion  current  (x)  collected  on  the  negative  electrode  of  radiois  f 


is  equal  to  all  of  those  ions  generated  in  the  cathode  sheath,  or 

.Pz  y  b(n) 

T  =  Z77-<S  /  r.Af,  f  . 


(35) 


o 
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Since  b(r2)  the  electrode  separation,  the  current  in  Eq  .  (35)  is  seen  to 

saturate  when  b(r2)=d  for  all  r^  • 

If  the  ion  generation  rate  is  unifom  over  the  sheath,  1  .e  .  S{r^,z^)=S, 
then  -rom  Eq.  (34)  V-  V^=(eS/[i+  6^) (b^/z ) ,b=  r/S=l/AeS  and  V=[4sV  e^e^j'^/^ 
(i/a)  which  is  identical  to  Eq.(29). 

Tlie  dependence  of  current  on  voltage  given  by  Eqs  .  (34)  and  (35)(or  by  Eq. 
(29))  will  break  dovm  when  diffusion  terms  become  significant.  In  addition, 
for  large  values  of  E(z)/p  in  the  sheath,  the  electrons  can  acquire  sufficient 
energy  to  cause  appreciable  multiplication  of  the  current  by  electron-atom 
‘^^llisions  and  this  effect  is  likely  to  be  of  particular  importance  in  a  gas 
mixture  such  as  neon-argoi. 

Computer  Programs  for  Ion  liansport  Model 

The  solutions  of  Eqs.  (34)  and  (35)  bo  obtain  the  current -voltage  charac¬ 
teristic  requires  numerical  as  well  as  iterative  techniques  .  Two  programs  were 
written  for  solution  on  the  IIM  7094  computer .  The  first  program  was  written 
for  a  simplified  single  integral  expression  for  b(r^)  and  the  program  details 
are  described  in  another  report .  The  second  program  was  written  for  the 

complete  Eqs.  (34)  and  (35)  and  solved  the  triple  integral  expression  Eq.  (34) 
for  b(r2).  Both  of  these  versions  are  described  briefly  and  a  comparison  of 
the  computed  I-V  curves  is  given. 

Simplified  Single  Integral  Version  In  Section  C  of  Reference  (l)  an 
integral  equation  for  b(r2)  for  the  cathode  sheath  thickness  was  witten  by 
analogy  with  the  uniform  ion  generation  rate  model .  For  a  circular  parallel 
plate  ionization  tube  of  radius  P2,cm  and  gap  d,cm  in  which  one  of  the  elec¬ 
trodes  contains  a  uranium  disc  of  radius  p^,cm,  the  integral  equation  in 
Reference  (l)  for  b(r2)  reduces  to  two  expressions  depending  upon  whether  the 
cathode  sheath  is  at  the  uranium  or  bare  electrode  .  Those  expressions  are 


b(^c)  ^ 

,  i'/3 

sd) 

(36a) 

i 

‘j 

and 

b?kj  ..  , 

1 

d  )  . 

4 

(36b) 

- 

.  Zt  =  c/  -  _ 
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The  tube  current  Is  ottalned  from  Eq.  (35)  using  bCr^)  from  Eq.  (36),  i.e. 


and 


Pu 

p^rre,  j  j  oi^ 

^z,-d-b(KO 


(37b) 


where  b  (r^)  gives  the  sheath  thickness  and  the  tube  current  when  the  uran¬ 
ium  electrode  is  negative  (cathode  sheath  at  the  uranium  electrode)  and  b"(r  ) 
gives  the  sheath  thickness  and  I  the  tube  current  when  the  bare  electrode  is 
negative  (cathode  sheath  at  the  bare  electrode).  Equations  (36)  and  (37)  were 
coded  in  Fortran  II  for  solution  on  the  IM  7094  computer.  All  numerical 
integrations  were  performed  by  Simpson's  rule,  using  a  fixed  number  of  steps 
(input)  for  each  variable.  Hie  solution  of  the  integral  Eq .  (36)  for  b(r,  ) 
required  iterative  techniques  and  this  was  done  with  an  existing  iterative 
Fortran  subroutine.  The  input  for  S(r2,Z2)  was  provided  as  a  tabular  function 
from  the  ion  generation  rate  code .  Values  of  the  current  l'^  and  I"  were  com¬ 
puted  for  a  given  voltage  range  at  the  voltage  points  fixed  by  the  input  val¬ 
ues  of  the  upper  and  lower  bounds  and  the  increment  size  .  The  computed  values 
of  the  current  were  printed  out  in  tabular  form  for  each  voltage  point  and 
these  data  could  also  be  punched  out  in  card  form.  Also  a  General  Motors 
Research  Laboratories  raultiplotting  subroutine  (which  uses  a  CRT)  was  pro¬ 
grammed  into  the  code  so  that  the  complete  I-V  c^urve  could  be  obtained  in 
graphical  form  by  itself  or  together  with  the  experimental  I-V  curve  if  that 
data  had  a] so  been  provided  as  card  input . 


^iple  Integral  Version  The  single  integral  expression  Eq.  (36)  for  b(r  ) 
gave  only  an  approximate  solution  to  Eq,  (3^) .  Tne  complete  solution  of  Eq.^ 
(3^)  foi"  ^(r^)  was  accomplished  numerically  by  a  combination  of  triple  inte¬ 
gration  and  numerical  techniques  .  Tliis  program  was  witten  as  a  subroutine 
and  then  substituted  into  the  main  control  program. 
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Discussion  of  Solutions  for  Current -Voltage  Characteristics 

Examples  of  the  code  solutions  for  the  current -voltage  characteristics 
will  be  given  in  the  next  section  on  the  analysis  of  the  inpile  data.  The 
discussion  here  will  be  limited  to  a  comparison  of  the  output  from  the  single 
and  triple  integral  versions  for  the  ion  transport  model  and  to  a  comparison 
of  the  I-V  characteristics  using  the  linear  and  quadratic  energy  loss  models. 

Comparison  of  Ion  Transport  Solutions  Code  solutions  for  the  I-V 
characteristics  were  obtained  using  both  the  single  and  triple  integral  solu¬ 
tions  to  Eq .  (34)  for  the  case  of  2h0  torr  xenon  in  a  neutron  flux  of  1  .txlO^^ 
cm  sec  .  The  maximum  deviation  over  the  entire  voltage  range  (-100  to  lOOv) 
was  less  than  2/o.  For  this  cast  the  sheath  thickness  is  very  small  (-^.03  cm) 
and  therefore  the  variation  in  S  across  the  sheath  is  very  small.  In  order  to 
make  the  comparison  with  a  large  sheath  thickness  (b(r2)  .3  cm)  runs  were 

made  on  each  version  of  the  code  for  the  case  of  240  torr  neon  at  a  neutron 
flux  oi  1 .2x10  cm  sec  .  The  maximum  deviation  in  this  case  was  less  than 
'^5/°*  These  data  showed  that  the  single  integral  equation  was  a  good  approxi¬ 
mation  for  the  sheath  thickness  b(r2'  .  The  triple  integral  solution  did  not 
have  an  excessive  running  time  and  so  it  has  been  used  throughout  the  data 
analysis  . 

Comparison  of  Energy  Loss  Models  Code  solutiono  for  the  I-V  character¬ 
istics  were  obtained  using  both  the  linear  and  quadratic  energy  loss  models 
for  the  case  of  240  torr  xenon  in  a  neutron  flux  of  1 .2xl0^^cm"^sec‘^ .  The 
curves  are  shown  in  Fig.  9*  As  pointed  out  in  Section  II  some  trouble  was 
experienced  with  the  numerical  solution  for  the  triple  integral  solution  for 

^2 — curves  in  Fig.  g,  therefore^  are  for 
a  negative  applied  voltage  where  the  sheath  is  at  the  electrode  opposite  the 
uranium  (z^-d)  and  the  quadratic  energy  loss  solution  can  be  obtained  without 
the  integration  problem.  The  current  values  for  the  quadratic  energy  loss 
model  were  about  9°/o  lower  than  the  values  for  the  linear  model  and  this  is 
consistent  with  reduced  ion  generation  rate  found  for  the  quadratic  energy 
loss  model  in  Section  II. 

Since  this  difference  in  current  values  for  the  two  energy  loss  models  is 
within  the  accuracy  by  which  some  of  the  other  variables  are  known^  such  as 
neutron  flux, the  simpler  and  more  rapid  single  Integral  solution  for  the  linear 
energy  loss  model  was  used  in  all  further  computations . 
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V  (volts) 


IV.  AIMLYSIS  OF  INPILE  DATA  ON  IONIZATION  TUBE 


The  inpile  data  obtained  ^vlth  the  ioni nation  tube  on  pure  and  mixed  noble 
gases  were  described  in  detail  in  a  previous  report. A  sketch  of  the  ioni¬ 


zation  tube  is  reproduced  in  Fig.  10. 
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EUCrRODCS 
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u  backup  Disc 
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HUIARC  k'.UD 


high  density 

AlUMINA 


Fig.  10.  Ceramic-metal  ionization  tube 
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The  tube  had  3  Kovar  electrodes,  essentially  planar,  which  were  brazed  to 
ceramic  rings  to  form  two  ionization  chambers  each  of  gap  O.3  cm.  In  the  ioni¬ 
zation  experiments  described  here  only  the  fission  fragment  ionization  process 
was  of  interest  so  thermionic  electron  emission  was  suppressed  by  using  a  thin 
uranium  foil  which  minimized  heat  generation.  Tlie  upper  chamber  did  not  con¬ 
tain  uranium  and  ionization  ^•ra.s  due  to  7  irradiation  only,  whereas  the  ioniza¬ 
tion  produced  in  the  lower  chamber  v/as  due  to  both  fission  fragments  and  7  radi 
ation.  By  making  current -voltage  measurements  with  both  diodes  it  was  possible 
to  establish  that  the  ion  generation  rate  due  to  7  radiation  ^■ras  only  l^o  of 
that  due  to  fission  fragments  . 

Current -voltage  characteristics  for  various  gas  pressures  of  neon,  argon, 
xenon  and  neon-argon  (Ar/Ne=10  were  obtained  at  reactor  power  levels  of  10, 
100,  1000  and  2000  kilowatts.  These  data  have  been  analyzed  with  the  ion  gen¬ 
eration  rate  code  using  the  linear  energy  loss  model  and  the  current-voltage 
characteristic  code  using  the  triple  integral  solution  for  the  sheath  thick¬ 
ness  .  The  values  of  the  constants  used  in  the  equations  are  listed  in  Table  2 
The  fit  of  theory  to  experiment  is  discussed  for  the  pure  gases  first  and  then 
for  the  neon-argon  mixture  . 
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TABLE  2 
Constants 


Average  u-roup  Properties  of  the  ■t’lGslon  Frag;ment s 


Initial  kinetic  energy,  E  Mev''"^ 

( 1  O 

Mass,  Amu'  ‘ 

Initial  charge,  Z^, 

Range  in  air,  R  , 

3.’ 

Range  in  U-Ni,  R^,  cm 

Energy  to  create  an  ion  pair,  u,  ev 
w/w^  where  w=(  ' 

w(Ne)  (v^=36.8  ev^lO)) 
w(Ar)  (w^=26.)! 
w(Xe)  (v/-  =21.9 

L4<i  ‘  \ 

w(Ne-Ar,0  .l/o  Ar'(w  =26.1'^^^) 


Light  Fragment 
(j=l) 


6.65x10’ 


1.06 

39.0 

28.0 

23.2 

27.7 


Heavy  Fragment 

(.1  2) 


5  .05x10’ 


1 .11 


40.8 

29.3 

24 .3 
29.0 


Fission  Cross  Section,  Q^,  cm’ 


21 .7 


Properties  of  the  Gases 

Ion  mobility  at  76O  torr  and  273  °K,  |i  ,cm^volt‘^sec'^ 
p^(Ne  )  in  Neon^^^^ 
p^(Ne^)  in  Neon 
p^(Ar  )  in  Argoi/^^^ 

^^(Ar^)  in  Argon 
p^(Xe  )  in  Xenon^'^^^ 

|j,^(Xe^)  in  Xenon 
p^(Ar^)  in  Neon^^^^’^"^ 
ti^(Ar2)  in  Neon^^^  ^ 


1.9 

.58 


Recombination  Coefficient  (Dissociative) 
a  (Ne^),  cm-^sec'^ 

,  .  .(16) 


0^r(^^'p.'  cm^sec"^ 


C^,  (Xep),  cm^sec  ^ 


,(16) 


2.2  X  10"  ' 

6.7  X  10'^ 


1.4  X  10’ 
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Pure  Gases 


The  experimental  data  for  the  pure  gases  neon,  argon  and  xenon  at  2t0  torr 
pressure  and  a  neutron  flux  of  IxlO^^cm  ^sec  ^  are  shovm  in  Fig  .11 .  For  dimity  these 
data  are  shovm  as  points  although  each  I-V  characteristic  was  swept  out  with 
an  x-y  recorder  .  The  experimental  characteristics  generally  all  followed  a 
dependence  on  voltage  and  dependence  on  neutron  flux  as  predicted  by  Eq.(29). 

The  computed  I-V  characteristics  for  neon,  argon  and  xenon  at  240  torr 
pressure  and  a  neutron  flux  of  1x10  ^cm  sec  ^  are  shown  as  the  solid  curves 
in  Fig.  11.  The  fit  of  the  theory  to  the  experimental  data  is  good.  For  the 
negative  quadrant  the  entire  area  of  the  Kovar  electrode  was  used  for  the  ion 
collecting  area  while  lor  the  positive  quadrant  the  ion  collecting  area  was 
limited  to  just  the  uranium  area.  This  was  justified  by  the  follov/ing  argu¬ 
ments  . 

When  the  applied  voltage  is  negative  (Kovar  electrode  negative  with  respect 
to  the  uranium)  the  cathode  sheath  is  at  the  Kovar  electrode  and  the  ion  gen¬ 
eration  rate  is  more  uniform  (see  Fig.  7)*  For  this  condition  the  fit  is  good 
and  the  theoretical  curves  for  neon,  argon  and  xenon  are  about  10 %  higher  than 
the  experimental  values  . 

For  a  positive  applied  voltage  the  cathode  sheath  is  at  the  uranium  elec¬ 
trode  .  The  ion  generation  rate  drops  sharply  with  distance  past  the  edge  of 
the  uranium  disc  (see  Fig.  7)  and  the  contribution  to  the  total  current  for 

rp^p  is  therefore  small.  More  importantly,  however,  for  r_>p^  the  Kovar 
.  cl 

supporting  the  uranium  bends  away  from  the  uranium  and  the  ions  must  be  trans¬ 
ported  to  the  electrode  acioss  a  region  where  the  ion  generation  rate  is  zero 
(see  -i' ig .  10)  .  It  can  be  shown  that  most  of  the  applied  voltage  occurs  across 

region  of  zero  ion  generation  rate  and  the  current  contribution  is  further 
reduced  .  The  code  was  not  programmed  to  take  into  account  this  detail  of  the 
tube  geometry  so  that  for  the  theoretical  curves  for  V  >  0,  the  ion  collecting 
electrode  radius  Pp  was  set  equal  to  p^  .  Even  with  the  reduced  collecting 
area  for  V^O  the  current  was  higher  with  the  cathode  sheath  at  the  uranium 
electrode.  This  results  from  the  decrease  of  the  ion  genera »-.ion  rate  with 
distance  from  the  uranium  surface  (see  Fig.  5). 

The  shape  of  all  of  the  experimental  I-V  curves  for  the  pure  gases  fol¬ 
lowed  very  closely  the  V^^^  dependence  predicted  in  Eq.  (29).  The  increasing 
current  values  (at  constant  voltage)  for  the  gases  in  the  order  neon,  argon 
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and  xenon  reflect  the  increase  of  fission  fragment  energy  loss  and  the  increase 
of  the  ion  generation  rat;e  with  higher  atomic  mass  of  the  gas  atom.  The  com¬ 
puted  ion  generation  rates  at  ihe  midpoint  of  the  tube  are  listed  at  the  bot¬ 
tom  of  Fig.  11.  For  neon,  argon  and  xenon  they  are  (l.8^^.6  and  8  .l)xlO^*^cm" 
sec  ,  respectively. 

Neon -Argon  Mixture  \^^,th  Ar/Me^l  x  10~^ 

The  experimental  Z-V  data  for  the  neon-argon  mixture  at  2k0  torr  and 
13  -2  -1 

1 .0x10  cm  sec  neutron  flux  are  plotted  as  the  points  in  Fig .  12 .  For 
!v|c25  volts  the  current  is  roughly  propoi'tional  to  (see  dashed  curve) 

and  is  about  tvrice  that  for  pure  neon  (Fig.  ll)  .  At  about  25  volts  the  elec¬ 
tric  field  becomes  sufficiently  gi-eat  to  cause  electron  impact  ionisation  in 
the  sheath  and  the  current  then  increases  rapidly  with  voltage  . 

The  theoretical  curve  ■was  obtained  with  the  code  using  w  =26.1  ev/ion 

■  (11)  ^ 

pair  for  the  neon -argon  mixture  viith  Ar/Ne=10  and  the  range  of  the 
fission  fragments  used  v/as  the  same  as  for  pure  neon .  The  experimental  cur¬ 
rent  was  higher  than  the  theoretical  current  for  (-25  volts <=¥*== +25  volts) 
and  this  was  attributed  to  the  diffusion  of  Ar  ions  from  the  bulk  of  the 
plasma  into  the  sheath.  Tliis  conclusion  is  supported  by  the  fact  that  good 
agreement  was  obtained  at  hOO  torr  pressure  (as  sho\m  in  Fig  .  I3  for  the 
negative  quadrant  when  the  cathode  sheath  is  at  the  Kovar  electrode)  and  as 
the  pressure  decreased  down  to  30  torr  the  observed  current  became  progres¬ 
sively  greater  than  tiiat  computed  .  The  ion  generation  rate  for  the  neon-argon 
mixture  at  2t0  torr  (2 .6xl0'^*^cm  ^sec  ^)is  t07o  greater  than  that  for  pure  neon 
and  this  is  attributed  to  the  Penning  effect  vMch  results  from  the  ioniza¬ 
tion  of  argon  by  the  metastable  states  of  neon. 

The  good  agreement  between  theory  and  experiment  for  pressures  of  2h0 
and  400  torr  is  taken  to  indicate  shat  the  computation  of  the  ion  generation 
rate,  including  the  geometrical  effects,  is  valid.  Table  3  gives  values  of 
S  computed  for  the  midpori.nt  along  t.he  tube  axis  for  the  gases  studied  .  These 
values  of  S  are  higher  than  the  average  values  previously  reported from  a 
one  parameter  fit  of  t,he  ion  transport  model  to  the  experimental  data.  The 
present  calculations  for  S  using  the  computer  code  can  now  be  utilized  to 
determine  the  spatial  distribution  of  S  for  a  irlde  range  of  parameters  . 


Also  shown  in  Table  3  values  of  electron  density  corresponding  to  the 
ion  generation  rates  .  The  electron  densities  for  the  pure  gases  were  calcu¬ 
lated  from^the  ion  generation  rate  and  the  dissociative  recombination  coef¬ 
ficient;  i  .e  .  N^=(s/q;^)' .  'The  determination  of  the  electron  density  for 
the  Ne-Ar  mixture  is  discussed  in  Section  C  of  this  report . 


Neon 


Argon 


Xenon 


TABLE  3 

Ion  Generation  Rate  and  Plasma  Density  at  Midpoint  of  Gap 
Gas  pressure  =  240  torr 


Neutron  flux 


Ion  Generation  Rate 

(  -3  -1\ 

_ ycm  sec  ) _ 

.  > 

1.8  X  10^^ 

3 .6  X  10^^  ^ 
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Fig.  11.  Current -voltage  data  in  neon^  argon  and  xenon. 


Fig.  12.  Current -voltage  data  in  the  neon-argon  mixture. 


Kma) 


0  5  10  15  20  25  30  35  40 

V  (volts) 


Fig.  13,  Current -voltage  data  in  Ne-Ar(0.iyo  Ar)  at  400  torr  pressure. 

The  applied  voltage  is  negative  and  the  cathode  sheath  is  at 
the  Kovar  electrode . 
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SECTION  B 


REACTION  KINETICS  OF  THE  FISSION  FRAGMENT  PLASMA* 


ABSTRACT 


atomic  ^  and  lom  c  reactions  in  a  plasma  formed  by  fission  fragments  in  Pen- 
ing  gas  mixtures  are  analyzed  .  Seven  equations  are  witten  for  atomic  and 
molecular  ions  of  both  major  and  minor  species,  metastable  atoms  of  the  major 

other  than  metastable  of  the  major  species,  plus  charge 

conS?tnS;  fl’Nf “f ®  solutions  are  obtained  for  the  major  and  minor  gas 
constituent  density,  to  give  maximum  electron  density.  The  classical  Penning 

+  source  rate  because  of  conversion  of  major  species 

IrrimTOrtanris  tl  f  d”  ^  relatively  small  advantage.  Much 

rate  IS  the  decrease  in  minor  species  atomic  ion  loss  rate.  This 

atomirions  0?^^^  operating  in  the  regime  of  parameters  where 

omic  10ns  of  the  minor  species  must  form  molecular  ions  via  3-body  collisions 
This  process  can  be  made  very  slow.  Examples  are  given  for  neon-argon  gas^ 
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1 .  The  Problem 


We  consider  a  tube  such  as  that  discussed  in  Section  A,  filled  \^^ith  either 
a  single  gas  species  or  a  Penning  mixture  and  containing  a  uranium  fission 
source  such  that  when  operated  in  a  reactor  a  copious  flux  of  fission  frag 
ments  is  passed  through  the  gas  .  The  fragments  may  ionize  or  excite  neu¬ 
tral  atoms;  metastable  atoms  may  be  formed.  These  particles  interact  among 
themselves ^ and  with  each  other ^  can  diffuse  to  the  tube  walls  and  the  ions 
can  recombine  with  electrons  (there  are  no  negative  ions).  Despite  the 
appreciable  variation  in  the  fission-fragment  ionization  rate  across  the 
tube  discussed  in  Section  A,  the  effect  of  such  variations  is  small  in  the 
present  study  of  the  number  densities  of  various  species  in  the  plasma. 

Hie  particle  density  distributions  in  the  plasma  are  dominated  by  the 
influence  of  such  loss  processes  as  diffusion  or  recombination  and  the  ion 
generation  rate  can  be  considered  effectively  uniform  throughout  the  tube 
volume  without  much  error . 

The  problem  is:  firsts  to  create  equations  governing  the  processes;  second^ 
to  find  coefficients  for  the  rate  processes;  thirds  to  find  solutions  for 
the  interacting  particle  densities. 


2 .  The  Equations 

Specifically^  we  assume  a  majority  gas  of  fixed  neutral  density  N  ^  which  is 
ionized  to  an  ion  density  .  There  can  also  be  metastable  atoms^vLth  den¬ 
sity  excited  states  that  are  not  metastable ^  and  even  molecular  ions 
N2+ .  There  may  be  a  minor  Penning  constituent  with  neutral  density  A  ^  atom! 
ion^  density  A_^^  molecular  ion  density  A^,^ .  We  ignore  the  presence  of^all 

excited  states  A^  and  A^^  because  these  cannot  yield  ions  via  any  linear  pro¬ 
cess  . 

The  first  equation  is  charge  neutrality 


"e  -  t  -  ”2.  -  A  -  =0  (1) 

Tlie  six  equations  for  N  ,  M  ,  ,  A_^  and  A  are  sho\\rn  on  the  next  page 

in  semi-tabular  form.  The  steady  state  solution  with  all  derivatives  zero 
is  desired  .  We  now  discuss  these  equations  one  at  a  time 
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Equations  for  Particle  Density  Time  Derivatives 


(a)  N_^ 


N_^  is  made  by  a  volume  source  route  S  N  ;  note  that  So  is  the  source  rate 
density  per  neutral  atom^  and  has  the°dimension  (sec"^)  .  Tliese  ions  diffuse 
(coefficient  D^)  and  can  recombine  with  electrons  (2-body  radiative  Ci, 

3-body  with  a  neutral  C2,  3-body  with  extra  electron  C^)  .  We  shall  neglect 
C3  because  Nq  »  hence  C2NQ  >-3>  .  Note  that  the  various  coefficients 

C  have  various  dimensions  suitable  to  each  occasion. 

N+  is  also  produced  by  collision  of  two  metastables  (C^)^  and  are  converted 
via  3“body  molecular  ion  fomiation  (Ci)  to  N24- •  They  are  also  converted  to 
via  2-body  charge  transfer  with  A^fcg) .  We  shall  ignore  the  3 -body  col¬ 
lisions  C^N+AqNq  compared  with  . 

(b)  N 

X 


N^  is  created  by  primary  radiation  (S^)  and  from  a  fraction  Co  ^ 1  of  the  N 
recombinations.  It  is  lost  via  radiation  i'^x)  >  by  the  molecular  ion  con¬ 
version  process  Co  (the  excited  atom  strikes  a  neutral  to  form  a  molecular 
ion^  before  radiating). 

(c) 


\  is _ created  by  primary  radiation,  by  a  fraction  C^q  <  1  of  the  N+  recom¬ 
binations,  and  by  a  fraction  C22_  <C  1  of  the  N2+  recombinations.  It  is  lost 
via  diffusion  Dm,  and  via  metastable -metastable  collisions  Cp2;  we  must  have 
^2.2  ^  Cy .  It  is  also  lost  via  2-body  destruction  with  neutrals  (Cpn)  and 
3"body  destruction  Cpt .  Very  importantly,  N^j^  is  converted  to  A^.  via  the 
Penning  process  . 

(d)  N2+ 

N2+^(and  A+  and  A2+)  is  nob  created  in  appreciable  amount  via  the  primary 
radiation.  Thus  N2+  is  created  by  the  3-body  conversion  of  N^  and  conver¬ 
sion  of  Nx,  all eady  described.  It  is  lost  via  diffusion  (D2),  dissociative 
recombination  Cpg  and  charge  transfer  with  Aq  to  form  A+  . 

(e)  A^ 


A^  is  formed  most  importantly  via  the  Penning  process  Cp,-,  and  also  via  Cr 
and  CpY^ already  described.  It  is  lost  via  diffusion  D^, ^radiative  recombin¬ 
ation  with  electrons  (Cp0),  3-body  recombination  using" neutral  N^  (Cpn),  and 
3-body  collisions  to  form  (C20)  •  Note  that  this  last  process  reads 


the  collisions 

A,  t  N  +  N  - ^  ? 

+  00 

occur  much  more  frequently,  but  we  deny  the  existence  of  (AN)  . 

W)  Ag, 


A2+  is  formed  via  C2Q  just  described,  and  is  lost  by  diffusion  (D.  )  or  dis¬ 
sociative  recombination  C21. 
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3  •  A  Pure  Gas  (Argon) 


The  preceding  eq^uaiions  are  too  numerous  for  comfort  ^  so  we  seek  preliminary 
insight  from  a  simple  case.  We  assume  only  one  gas,  thus  eliminating  Eqs . 

6  and  7*  The  principal  effect  of  was  to  make  A^,  which  now  can  no  longer 
happen.  The  only  remaining  contribution  could  be  N4.  via  and  this  is 

non-linear  and  (if  Sq  ^  Sq)  cannot  lead  to  any  large  change  in  .  Thus  we 
ignore  Njjj  also.  Finally,  we  shall  ignore  on  the  following  more  subtle 
basis:  If  Nq  is  high  enough  to  make  an  important  process,  then  Nq  is 
probably  also  high  enough  to  make  a  very  important  process  too .  Thus 
virtually  all  the  ions  are  molecular;  but  surely  Sq  1>  S,  ,  and  many  states 
radiate  in  spite  of  this  process .  Hence  the  contribution  of  Cn  is  not 
expected  to  be  large .  ^ 

Finally,  we  replace  V  by  l/ ,  where  A  is  the  diffusion  length.  This 
approximation  is  valid  in  the  regime  where  diffusion  dominates,  i.e.,  where 
the  term  is  important .  It  is  incorrect  only  in  the  regime  where  it  is  small 
compared  to  some  other  loss  term.  The  resulting  steady  state  equations  are: 

0  =  ,  (9) 

A" 

0  =  .  (10) 


We  are  interested  in  the  density  N^,  and  the  effect  of  varying  the  param¬ 
eters  .  Particularly,  we  are  interested  in  the  effect  of  varying  N  .  Now 
D  a  1/Nq,  hence  we  write  ° 


”1  = 
'>2  = 


\ 


(11) 


and  appears  explicitly. 

There  are  two  simple  limits  for  these  equations.  At  low  pressure  N  - >  0, 

diffusion  is  large,  3-body  conversion  is  negligible,  and  all  ions  are  N4. . 

The  solution  is 


N 

e 


S  N 
00  ' 


The  other  limit  obtains  at  very  high  pressure  -  no  diffusion, 
sion  to  N^^,  and  a  density  limit  given  by  recombination: 


(12) 

instant  conver- 


(13) 
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nder  some  circumstances^  there  is  an  interesting  anomaly  in  the  mid-range 
describable  physically  as  follows  *  If  the  recombination  Cpg  is  large 
(circumstances  can  always  be  arranged  so  that  it  is  "large"),  we  wish  to 
prevent  ions  from  recombining.  But  if  the  Ng  is  large,  this  can  only  be 

done  by  preventing  formation  of  Np^  in  the  first  place .  Thus  we  can  set 
Ng  =  and 


N 

e 


S  N 
o  o 


''l  2 

■2:7 " 


at*, 


c 


This  quantity  is  a  maximum  at  the  neutral  density 


(lU) 


where 


0  1/3 

(2  V  /\\) 

(15) 

p  0  1/  3 

Sjh  /ay 

(16) 

It  is  instructive^  to  compare  Eqs .  I3  and  I6 .  Since  Ne  a  N  at 

principle  make  Ng  arbitrarily  large  just  by  increas 
ing  Nq.  But  this  IS  unattractive,  because  of  plasma  resistance,  and  so 
lorth.  If  Nq  is^limited  by  practical  considerations,  then  Eqs.  I5  and  I6 

may  in  fact  predict  a  more  interesting  regime  of  operation  at  relatively 
low  density  ^ 


Of  course,  there  is  nothing  wr  ng  with  solving  Eqs .  8  -  10  directly.  We  have 


N  ^  +  N 
e  e 


We 


K, 


S  N 
o  o 


N 

0 


S  N 
o  o 


K, 


2 

+  C,  N 
4  o 


'1 

^  b  +  C,  N  ^ 

2  4  o  i 


^1  2! 
- o  +  0,  N 

N  ^  o 

o' 


0  . 


(17) 


but  the  behavior  of  tlie  solution  is  not  so  readily  apparent  without  the 
physical  argument . 


We  illustrate  the  situation,  using  argon,  and  take  the  following  coefficients: 

^0^0  (1*65)  10^^  (cm"^  sec“^)  at  240  Torr, 

so  that 


=  1.96  X  10" 3  sec~^ 

^16  =6.5  X  10’^  cra^  sec“^. 
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For  the  coefficients  K,  we  note  a  mobility  1 .6  for  the  atomic  ion  (here  called 
N_^) .  We  use 


2  k  T^' 

(  "0) 

\  e  - 

e 

N 

0 

and  set  kT^/e  =  0.08  ev.  Since  1-i^.N  =  constant  and  =  lo6  at  N  = 
2.79  X  I0I9  cm“^,  we  find  ° 


=  0.71  X  10^^ 

19  ? 

Similarly  K2  =  1.15  x  10  for  the  molecular  ion,  using  |i  =  2.6  cm^/volt-sec . 

Finally,  we  must  guess  about  .  It  is  reported  that  the  frequency  for 

He'*'  +  2He — >He2  is  mi 'JO  p2^  where  p  =  Torr .  We  shall  take  100  p2,  whence 
we  find 


C|^  =  8  X  10~^^. 

Note  that  if  C|^  is  smaller,  the  anomalous  electron  density  will  be  higher. 

Figure  1  shows  the  result  of  solving  Eq.  I7,  for  three  flux.  One  is 
So  =  1.96  X  10"3  sec-1,  corresponding  roughly  to  10l3  neutrons/cm^sec,  and 
the  other  two  are  flux  of  0.1  and  10.0  times  that  flux.  The  previous  popu¬ 
lar  operating  density  was  a/ 8.5  x  10I8  cm-3,  -  2U0  Torr.  The  region  around 
2.8  X  10  -T  cm-3,  8  Torr,  looks  interesting.  In  particular,  the  anomaly  is 
large  at  high  flux.  It  appears  that  the  simple  eqs .  I5  and  I6  are  valuable, 
because  they  predict  the  maxima  shown  by  crosses  in  the  three  cases .  The 
approximation  is  good  when  a  maximum  exists,  as  we  would  expect. 

•  Simplified  Description  of  the  Penning  Mixture  Physics 

Figure  2  shows  all  the  reactions  given  in  Eqs .  2-7  •  The  electron  population 
is  not  shown  specifically,  but  the  appearance  of  electrons  or  atoms  with  any 
reaction  is  shown  by  an  appropriate  input  arrow.  Note  that  the  only  feed¬ 
back  loop  is  CY(meta stable -met a stable  collisions  leading  to  ionization)  . 

These  equations  can  easily  be  programmed  for  a  computer  but  here  we  desire  a 
little  preliminary  insight . 

Figure  3  shows  in  bold  outline  two  parts  of  the  reaction  paths .  The  upper  part 
is  the  pure  gas  problem  described  in  the  previous  section.  We  turn  our  atten¬ 
tion  now  to  the  whole  diagram;  then  in  a  few  moments  to  the  lower  emphasized 
part . 

The  object  of  having  a  Penning  mixture  is  to  increase  the  charge  density  in 
the  plasma.  We  can  do  this  by  increasing  the  ion  production  or  decreasing 
the  ion  loss.  The  Penning  mixture  increases  the  ion  production  by  converting 
metastables  of  the  majority  gas  to  ions  (of  the  minority  gas).  Thus  in  Pig .2 
(or  3),  we  add  the  source  S2  to  the  source  Sq  ,  But  ^2  ^  ®o>  ^hile  we 
increase  the  ion  production,  we  do  so  only  modestly. 

Note  that  the  huge  Increase  in  ionization  normally  thought  to  be  associated 
with  a  Penning  mixture  is  here  absent.  The  reason  is  simple:  the  classic 
benefits  of  a  Penning  mixture  come  about  from  having  an  electron  energy 
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Figure  2  -  Reaction  diagram  corresponding  to  Eqs  .2-7,  for  production 
and  loss  of  ions  N+,  N2^,  A^,  and  A^^ . 
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Figure  3  -  Principal  reaction  sub-diagrams  for  ions  (bold  dotted) 
and  for  ions  (bold  dashed). 
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distribution  f(u)  of  the  shape  shown  by  curve  1  of  Fig. 
applied  electric  field.  Because  of  "elastic  collision" 


4,  maintained  by  an 
energy  loss  at  low 


Energy  u 

Figure  4  -  ^tist's  conception  of  electron  distribution  functions 

(1)  in  which  electrons  are  heated  from  low  energy^  and 

(2)  in  which  electrons  slow  down  from  high  energy. 


energy,  some  electrons  have  energy  above  u^  and  make  metastables.  These  con¬ 
vert  to  ions  of  the  minor  species,  and  thus  we  have  a  good  supply  of  ions 
where  there  were  formerly  very  few.  The  electron  distribution  here  is  shown 
in  curve  2:  all  the  thermal  electrons  near  u  =  0,  plus  the  slowing  down  dis¬ 
tribution  at  higher  energy.  Here  ionization  proceeds  with  vigor,  and  the 
excitation  of  metastables  is  a  modest  correction.  In  summary,  the  classic 
Penning  effect  is  important  when  the  (energy/ion)  is  very  high.  In  our  cases 

the  ^ev/ion;  is  purposely  made  low,  and  only  moderate  Improvements  are  to  be 
expected . 

To  be  sure,  these  moderate  source  improvements  are  desirable,  but  much  more 
important  is  the  reduction  in  loss  rate,  a  quite  different  phenomenon.  In 
Figs.  2  or  3,  we  wish  to  block  the  paths  by  which  any  one  of  the  ion  species 
disappears,  without  unduly  hindering  its  production.  Indeed,  this  is  what 
happened  in  the  previous  section:  the  loss  rate  of  N.  was  minimized  to 
obtain  the  anomalously  high  electron  density. 

Now  look  at  the  lower  emphasized  part  of  Fig.  3.  The  production  terms  are 
different,  involving  the  metastable  but  here  we  see  that  the  loss  rate 
of  A+  is  C20A+A0N0.  Here  is  the  chance  to  reduce  the  ion  loss  rate  dras¬ 
tically.  Because  Aq<K  Nq,  an  effect  not  possible  in  a  pure  gas  is  at  hand. 

We  see  that  some  N4.  and  N2+  ions  can  contribute  to  A4.,  but  for  the  moment 
we  ignore  all  but  the  lower  emphasized  reactions.  In  justification,  we  point 
out  that  both  the  production  via  metastables  and  the  loss  of  A+  will  be 
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sensitive  to  the  ratio  Aq/Nq,  but  all  t>e  rest  of  the  variations  are  insensi- 
tive  to  this  ratio.  Thus  our  plan  is  to  explore  this  critical  production  loss 
chain  critically.  If  the  A+  density  is  desirably  high,  we  can  add  the  effects 
01  h+^and  ^2-1.  later,  because  they  will  be  relatively  small.  If  the  A+  den¬ 
sity  IS  disappointingly  low,  we  dismiss  the  whole  sequence. 


5 •  Simplified  Penning  Mixture  Equations 


The  discussion  of  Section  4  is  now  turned  into  arithmetic.  Equations  4,  6,  7 
and  1  become  respectively  ^  ^ 

2 


0  = 


K  N 

- p  -  C,  N  Aq  -  C,  N  -  C,  ,  N  N  , 

15  m  o  13  m  o  l4  m  o  ^ 


“nA' 


(18) 


0  = 


^15^0  "  rj  yy2  ■  ^ 


(19) 


0  = 


CppA  A  N 
20  +  o  o 


4  2+ 


N  A‘ 
o  ' ' 


2  ^ 


(20) 


0  =  N 


A  -  A 
+  2+ 


(21) 


Srn"'!  Sr  ion  recombination  C  o  and  C._;  and  we  have 

Equation  l8°yiSdr  described  in  the  sentences  ]ist  before  Eq.i; 


N  = 
m 


SqN 
2  o 


K 


(22) 


”  -2  "  =15^  *  =13“o  " 


NoA 


trivially,  hence  from  Eq .  19 


C,_S„(A  N  )N 
15  2'  o  o  o 


[5^2^  <=20(Vo)j  [-^2  ^ 


(23) 


In  obtaining  Eq .  23,  we  have  multiplied  above  and  below  by  N^,  in  order  to 

g  mSmL?!  stratagem  will  be 

^e  remaining  equations  can  in  fact  now  be  solved,  yielding  a  quadratic  form 
lor  1^0  .  It  IS 


K, 


No  +  N 


-  A, 


SA/r 


K, 


^oA 


(24) 


where  we  must  insert  A+  from  Eq .  23 . 
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While  suitable  for  an  elementary  computer,  Eq.  2k  will  not  be  very  tractable 
analytically  when  we  consider  maximizing  Ng  with  respect  to  both  Nq  and  Aq  . 
Thus  we  use  the  stratagem  that  led  to  Eqs .  l4  -  l6:  if  atomic  ions  are  effec¬ 
tively  blocked  from  converting  to  molecular  ions,  then  Eq.  23  will  represent 
Ng  to  a  good  approximation. 


Our  plan  is  to  maximize  Eq.  23  with  respect  to  N  and  Aq  as  well  as  possible. 
To  start  the  process,  suppose  for  a  moment  that  C13  and  C]_4  were  zero  -  that 
is,  there  is  no  metastable  destruction  by  2-body  or  3-body  processes;  i.e., 
all  diffuse  or  convert  to  A4. .  Equation  23  could  then  be  written  as 


where 


N  = 


C15S2X 


^3  ^20^ 


K  +  C.  c-y 
m  ly 


X  =  AqNq^ 


y  = 


(25) 


(26) 


The  variables  x  and  y  are  just  as  suitable  for  maximizing  N  as  were  A  and  N  . 

Now  we  see  from  the  form  of  Eq.  25  that  there  is  no  maximum!  In  fact,  N  - 

X— ^00,  y  “^0,  or  as  Nq-^oo,  and  in  a  suitable  manner.  Note, 

however,  that  the  physics  prevents  our  setting  A^  identically  zero.  This  result 
is  very  different  from  the  pure  gas  results  of  Section  3;  the  reason  for  this 
difference  is  important  to  us.  Here,  the  ionic  conversion  rate  is  proportional 
to  AqNq,  not  Nq'^  in  this  approximation.  Thus  we  are  at  liberty  to  adjust  Nq 
anyway  we  please,  keeping  AqNq  constant.  In  particular,  we  can  increase  Nq 
indefinitely,  thus  producing  additional  ionization  S2NQ,  and  reducing  diffusion 
to  a  negligible  value . 


It  is  evident  that  metastable  destruction  at  high  pressure  upsets  this  picture, 
and  we  must  reintroduce  either  C^o  or  *^l4^  both,  in  order  to  recover  reason¬ 
able  physics.  This  single  insertion  causes  formidable  complications^  because 
no  separation  into  two  separate  maximizations  can  be  found  (as  was  done  in  Eq. 
25).  The  best  method  of  treating  Eq.  23  exactly  seems  to  be  to  maximize  with 
respect  to  AqNq,  and  Nq.  Setting  bK/  ^  yields 


(AN  )‘ 
'00 


15  "20 


+  Ct_N 

13  o 


^  +  C  .  N  ^ 
14  o 


(27) 


Note  that  Eq.  27,  which  gives  the  mixture  ratio,  predicts  that  this  ratio  will 
be  a  strong  function  of  the  density  itself.  Setting  3a  /^N.  =  0  is  reported 
to  yield  ° 


A  N 
o  0 


3  K^K 

- ^  +  KC  N^  +  KC  n3 

yy^2  3  13  o  3  14  o 


''i3'^20A\^  "  2"i4''20 


3  15 


m  20  o 


(28) 


Equation  28  does  not  seem  useful  to  us,  and  has  not  even  been  checked,  in 
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principle,  the  ideal  Nq  is  determined  by  solving  Eqs .  27  and  28,  but  the  equa¬ 
tion  vail  be  nth  order  (if  we  have  =  0,  it  is  8th  order) . 

We  must  use  more  physical  reasoning  to  simplify  the  maximization  problem.  It 
^11  turn  out  that  the  metastable  diffusion  rate  is  small  compared  to  metastable 
destruction.  The  reason  is  that  in  Eq.  23  we  will  end  up  by  being  close  to  a 
balance^of  diffusion  and  A4.  3-body  conversion  (K^  vs  Coq)  .  But  the  3-body 
conversion^ is  very  unlikely  compared  to  any  2-body  collision  process.  Then  on 
the  supposition  that  K3  and  are  commensurate,  and  C15  or  Cn  are  reasonable 
-body  numbers,  we  must  conclude  that  metastable  diffusion  will  be  small  com¬ 
pared  to  other  metastable  disappearance  terms . 

(^e  more  fortuitous  circumstance  allows  us  to  complete  the  calculation:  Cni.^  0 

for  the  neon  metastables  in  neon  (see  Section  D)  .  With  these  simplifi¬ 
cations,  Eq.  27  yields 


N 

o 


"3^13 


°15^2oAV 


(29) 


In  addition,  Eq.  23  becomes  remarkably  simplified  by  inserting  Eq.  29: 


A,  = 


A^c  ^ 
^  ^5  ^20 


a3  +  k  n  2  12 

o  o  13 


A^c  ■* 
' '  15  ^20 


This  is  easily  maximized  to  yield 


Thus  finally 


■^3^13'  1'^^ 


N 


1*  -11/3 

A"Ci3C20J 


2C. _ 
o  15 


and 


A^(«WN)  =  ^^2 


2  .2 


15  A 


1/3 


(30) 


(31) 


(32) 


(33) 


(3^+) 
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As  an  example,  let  us  take  the  follovdng  values.  (See  Section  D  for  a  full 
discussion  of  the  appropriate  values  to  use.) 

jy  7-2  X  10^®  cm‘^  sec“^ 

=  3  *6  X  10  cm^  sec”^ 

=  1.8  X  10”^^  cm"  sec"^ 

^20  ~  25  X  lO”^^  cm^  sec’^ 

Sg  =  1.35  X  10"3  sec"^ 

We  take  f\=  0.1  cm,  and  find 


A  , 

o  -4 

=  1  X  10 

o 

s, 

Nq  =  3-9  X  10^®  cm"^  (a?  no  Torr) 

=  3  -1  X  10^^  cm"^ 

6 .  Effect  of  Other  Processes 

The  heavy  dashed  se<iuence  in  Fig.  3  represents  only  part  of  the  story.  We  now 
put  in  the  remainder  of  the  important  effects. 

(s-)  production  and  density. 

For  N^,  we  shall  use 

=  2.6  X  10”^ 

=  8  X  10"3^ 

=  5*52  X  10^®  cm"^  sec"^ 

To  first  order.  Eq.  lU  should  apply,  whence  at  =  3.9  x  10^®  cm"^,  we  find 
N+  8  X  10^  cm”^.  This  is  a  small  perturbation  to  the  various  densities. 
Nearly  all  the  ions  disappear  by  converting  to  N2+  via  process  C]^. 

(b)  — >A_^  charge  transfer  (C^) 

We  might  have  significant  additional  source  of  A  via  process  .  We  wish 
to  compute  the  ratio  6 


C.N  A 
6  +  0 

C_N  A 
15  m  o 


‘^15^ 


For  these  approximate  conditions,  we  have 


Using  Cg  =  3.9  X  10  cm  sec"  ,  we  find  this  ratio  to  be  10"-\  Thus 
2-body  charge  transfer  N4.  >  A+  is  not  important  in  this  instance.  Note,  how- 

calculation  is  made  here  for  specific  conditions  and  the  factor 
1  will  not  always  apply.  At  lower  pressure,  fewer  N4.  ions  disappear  by 
process  hence  more  will  disappear  by  charge  transfer  (if  diffusion  does 
not  take  over  meanwhile).  On  the  other  hand,  the  actual  ratio  is  so  small 

{10  0;  that  we  will  believe  pro  tern  that  process  will  never  be  important  in 
neon -argon .  ° 

(c)  density. 

Here  we  can  compute  Np4.  as  a  balance  between  production  S  N  ^  Ci  N 
and  recombination  CigHsf  .  For  neon  molecular  Ion  dlssocla?l?e  reoSmSina- 
tlon,  C15  =  2.2  X  10-f  cm3  sec-1,  and  we  take  N  =  3.1  x  10l2  from  the  previous 
calculation  of  .  Therefore 


for  the  density 
(d)  density. 


M  S  N 

^2+  ^  ^  ^ 


C-.  N 
16  e 


3.9  X  10^®  cm"3 


^  10^°  cm" 3 

,  and  this  is  negligible. 


We  need  not  even  calculate .  The  rate  of  formation  of  A^  and  N 
parable,  and  Thus  A^^  is  negligible. 

(e)  Recombination  of  A_^ . 

We  take  the  follo\4ng  numbers  for  C^g  and  (see  Section  D) . 

^18  4  X  10  cm^  sec"^ 


are  com- 


=  2.7  X  10"^^ 

We  wish  to  find  the  ratio 


^°18  ^  °19”o)  Ve 


Note  that  at  =  3.9  x  10^®,  the  2-body  radiative  and  3-body  Thompson  recom- 
bination  rates  are  about  equal.  Tnus  by  chance  we  are  in  a  reasonable  pressure 
regime  for  non-excessive  recombination  loss.  For  a  density  N .  =  N  =  3  1  x  10^^ 
we  find  this  ratio  0.1.  Thus  this  process  is  small  also.  ®  ' 

(f)  N  effects. 

n  N  contributes  only  to  Np  ,  an  already  negligible  density,  we  neglect 

all  effects  of  N  .  « ^  o 

2  ^ 

(g)  destruction. 

We  find  sj  10^®  -  10^^  cm  For  a  metastable  atom  to  find  another 
metastable  before  finding  an  argon  atom  would  require  the  cross  section  cor¬ 
responding  to  Ci2  to  be  ^  10-12  ^  V^e  have  no  sure  data  on  this  number 

at  tne  moment,  but  believe  that  it  cannot  be  that  large.  Thus  we  dismiss  the 
process  indicated  by  . 
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CONCLUSION 


7. 


Physically,  we  have  maximized  the  electron  density  N^  in  a  Penning  gas.  At 
these  conditions  on  majority  gas  density  Nq  and  minority  gas  density  A-  that 
give  maximum  Ng,  we  find: 

(a)  Metastables  Nj^  are  destroyed  by  conversion  to  A’^  and  by  non-prof i table 
destruction  in  the  gas,  in  about  equal  proportions. 

(b)  The  ions  are  almost  entirely  A^,  which  are  lost  via  diffusion  and  3 -body 
conversion  to  A24.  (which  disappear  forthwith);  these  two  processes  occur 
with  roughly  comparable  frequency. 

(c)  The  ratio  (N^  conversion  to  A+)/(Njjj  useless  destruction)  is  exactly 
equal  to  the  ratio  (A^  diffusion  1oss)/(A^  3-tiody  conversion  loss)  at 
the  optimum  conditions.  This  fact  can  be  seen  by  inserting  Eqs .  3I 
and  32  into  Eq.  23. 

Note  that  the  recombination  terms  C^g  and  C]_^  are  presently  negligible 
(  ^  0.1)  compared  with  other  A4.  loss  terms,  and  Ng  is  proportional  to  S2. 

^is  will  not  continue  indefinitely.  At  a  flux  about  10  times  higher 
lOlA  neutrons/cm"^sec)  the  recombination  terms  will  be  important,  and 
Ng  will  tend  toward  a  dependence.  The  corresponding  density  will  sub¬ 
stantially  exceed  10^3  cm" 3  however. 
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SECHON  C 


COMPUTER  SOLUTIONS  TO  THE  REACTOR  KINETICS  EQUATIONS 


ABSTRACT 


The  seven  simultaneous  equations  which  describe  the  reaction  kinetics  in  a 
plasma  generated  by  fission  fragment  ionization  of  a  Penning-type  gas  mixture 
were ^ solved  on  the  IB4  digital  computer.  The  results  of  the  first  parametric 
studies  on  the  effects  of  pressure^  concentration  of  the  minority  gas  species 
and  the  electrode  spacing  are  discussed.  The  first  runs  confirmed  the  approx¬ 
imate  solutions  obtained  analytically  in  the  previous  Section  B.  Using  the  full 
set  of  equations  in  the  code  it  was  found  that  in  the  recombination  of  the 
molecular  neon  ions  if  an  appreciable  fraction  (  .5)  produces  the  metastable 

states ^ of  neon,  these  in  turn  increase  the  production  rate  of  argon  ions  and 
significantly  increase  the  electron  number  density  ( x  2).  Reducing  the 
electrode  spacing  from  .33  to  .05I  cm  increased  the  diffusion  losses  and 
greatly  decreased  the  electron  number  density.  This  result  showed  that 
further  theoretical  work  must  be  done  before  this  tool  can  be  used  to  opti¬ 
mize  the  perfoimiance  of  a  thennionic  energy  converter. 
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IN'IRODUCTION 


Tlie  equations  governing  the  reaction  kinetics  of  the  fission  fragment  gener¬ 
ated  plasma  have  been  discussed  in  detail  by  Professor  D.  J.  Rose  in  the  pre- 
vious  Section’  p.  The  equations  are  summarized  in  Table  1  and  the  reaction 
diagram  is  reproduced  in  Fig.  1.  In  that  study  a  solution  was  found  anal^ci- 
cally  for  the  case  of  pure  argon  at  the  inteimiediate  pressure  where  a  maximum 
occurs  in  the  electron  number  density  versus  pressure  (  8  torr)  .  A  second 

solution  v/as  found  analytically  for  the  Penning  mixture  of  neon-argon  again 
for  the  case  where  a  maximm  occurs  for  the  electron  number  density  versus 
pressure  (  110  torr)  and  also  at  the  optimum  mixture  of  neon-argon  (A^/N  = 

10  ) .  In  both  of  these  solutions  a  number  of  simplifying  assumptions  were 

made  to  make  the  equations  tractable .  To  obtain  the  second  solution  only  the 
production  of  argon  ions  from  the  neon  metastables  which  were  produced 
directly  from  the  fission  fragments  (^2^0^  were  considered  .  The  presence  of 
the  neon  ions  was  neglected  as  well  as  the  production  of  neon  meta¬ 

stable  states  from  the  recombination  of  those  neon  ions  C  N  ).  The 

'  ll  16  2+  e 

effect  of  including  this  latter  term  was  not  discussed  but  it  would  increase 
the  source  of  ions  and  could  significantly  increase  the  electron  number 
density  for  0 .1  ^  —  1 .0 .  This  possibility  and  the  future  use  for  the 

A-Cs  system  pointed  up  the  need  for  a  solution  to  the  complete  set  of  equa¬ 
tions  v^nich  could  only  be  done  economically  on  the  I0i  computer.  Conse¬ 
quently^  a  computer  code  was  witnen  and  this  report  summarizes  the  i^esults 
of  the  initial  studies  with  this  code  . 


and  loss  of  Ions  N+,  A^,  and 
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Computer  Code 


A  problem  statement  was  submitted  to  the  GM  Computer  Technology  Department  for 
the  simultaneous  solution  of  the  seven  non-linear  algebraic  equations  .  Tlie 
differential  operator^  y  ,  was  replaced  by  1/a^,  whereby  is  the  diffusion 
length  and  the  diffusion  coefficient,  D^,  \,?as  replaced  by 

There  was  no  FORTRAN-II  subroutine  available  to  operate  in  the  GM  Research 
Laboratories  P-system  but  a  new  FORTRAN-IV  subroutine  developed  at  Purdue  was 
located  and  this  subroutine  was  operated  on  the  new  IM  "IBSYS"  system.  ¥. 
Butterfield  of  the  Computer  Technology  Department  vsrote  a  short  control  program 
for  this  subroutine  and  the  program  was  successful  on  the  first  test  problem. 

Conversion  Coefficients  for  Rate  Processes 

The  values  of  the  conversion  coefficients  used  in  these  first  parameter 
studies  with  the  "Reaction  Rate  Code"  are  listed  in  Table  2.  Some  of  these 
values  were  obtained  from  the  literature  and  some  were  computed.  They  are 
discussed  in  detail  in  Section  D. 

The  source  term  coefficient  Sq(N+)  was  obtained  as  outlined  in  Section  A  using 
the  fission  fragment  ionipfcion  code  for  neon  at  240  torr  and  at  a  neutron  flux 
of  <ti  =  1  X  10  cm  sec  (P  =  2  MW)  .  The  maximum  value  of  S  at  the  surface 

of  the  uranium  was  used  .  The  source  term  coefficient  for  the  neon  metastable 
production,  S2(Nj^),  was  obtained  from  the  above  value  and  tho  similar  source 
rate  value  from  a  neon-argon  code  for  the  fission  fragment  ionization  of  neon: 
argon  (l000:l)  at  the  same  pressure  and  at  the  same  flux. 

The  complete  equation  for  N^,  the  excited  states  of  the  majority  gas  species, 
is  built  into  the  code  but  the  source  term  S^  was  not  readily  detemined .  S^^ 
and  also  N^,  v/as  set  equal  to  zero  in  the  present  solutions  .  The  justifica¬ 
tion  for  ignoring  N  is  outlined  in  Section  D. 

X 

Parametric  Studies  with  Code 


The  parametric  studies  reported  here  are  not  complete  but  rather  vrere  meant  to 
check  out  the  code  and  also  the  values  calculated  in  Section  B.  Many  of  the 
simplifications  made  in  Section  B  are  also  contained  here  but  only  because 
some  of  the  coefficients  are  presently  not  knora .  Only  in  the  case  of  the  pro¬ 
duction  of  neon  metastables  from  the  'ecombination  of  the  neon  molecular  ions 
'^ll^q5^2+^e ^  was  it  thought  that  the  addition  of  a  neglected  term  could  appre¬ 
ciably  affect  the  electron  density. 

A  total  of  93  runs  (solutions  to  the  equations)  have  been  run  to  date.  The 
results  are  presented  in  graphical  form  in  the  follovang  sections  and  the  com¬ 
plete  solutions  to  six  of  the  runs  are  given  in  Table  3. 
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TABLE  2 


Conversion  Coefficients  for  Ne;A  Studies* 
(N  =  Neon  A  =  Argon) 


Source  Term  Coefficients 

=  2.6  X  10"^  cm"^ 
S,  (N  )  =0 

1'  X 

=  1-35  X  10"3  cm’3 

Reaction  Coefficients 


sec  ^  (unit  density) 

-1 

sec  (unit  density) 


=  (Varied)  0.0,  O.5,  1.0 


=  8  X  10  cm  sec 


=  3*9  X  10  cm^  sec"^ 


=  3 *8  X  10  cm^  sec"^ 


=  1.8  X  10  cm'^  sec^^ 
=  2.2  X  10  cm^  sec"^ 

=  3*9  X  10"^^  cm^  sec”^ 

I,  T  ^-12  3  -1 

=  4  X  10  cm'^  sec 

=  2.4  X  10  cm^  sec~^ 

=  2.3  X  10”^^  cm^  sec”^ 

=  6.5  X  10"^  cm^  sec"^ 


Diffusion  Coefficients 

V  I—  r“ 


=  5*52  X  10^  cm"^  sec"^ 

=  5 *5  X  10^®  cm"^  sec”^ 

=  8.74  X  10^®  cm”^  sec"^ 

=  7 *2  X  10^®  cm”^  sec"^ 

Kj^  =  6.8  X  10^®  cm”^  sec”^ 

Geometry  Coefficients 

A  =  d/jt  where  d  =  electrode  spacing 
A  =  0.1  for  most  of  these  studies 

Average  Lifetime  of  Excited  States 

Tx  =  1  X  10"  sec  (note  that  all  source  terms  of  N  =  O) 


more  detailed  evaluation  of  the  conversion  coefficients  is  to  be  found  in 
Section  D.^  Those ^ that  were  not  equal  to  zero  were  not  known  at  the  time  of 
writing  this  Section.  Hov/ever  the  influence  of  finite  values  for  th^m  has 
been  checked  and  found  to  be  insignificant  in  this  calculation. 
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TABLE  3 

Selected  Code  Solutions  to  Kinetic  Equations 

(Near  Maxima  in  Curves  of  N  vs  N  ) 

e  0 


Run 

No. 

X 

N 

10-ie 

0 

0 

< 

A 

N 

®-12 

X  10 

X  10-1° 

0 

1 — 1 

+ 1 

OJ  0 

M 

N 

""-12 
X  10 

X  10"^^ 

X  10 

10 .1 

3-9 

1 

0 

1 — 1 

X 

0.0 

0.1 

2.918 

.8331 

1.578 

.2483 

2.893 

5.800 

11.3 

3.9 

1 

X  10“ 

I, 

0.5 

0.1 

5.393 

.8331 

.8542 

.4873 

:;.0(5 

5.831 

14 .4 

4.0 

1 

X  10"^ 

1.0, 

0.1 

7.664 

.8123 

.6165 

.7266 

7.849 

6.142 

13.2 

2.0 

1 

X  10"^ 

X  10“^ 

\. 

0.5 

0.1 

3.192 

1.620 

.7367 

.1216 

3.167 

15.26 

17 .7 

8.0 

1 

0.5 

0.1 

3.154 

.4062 

2.997 

.6994 

3.120 

2.435 

19.9 

14.0  1x10  0.5 

Effect  of  Pressure 

.0162 

1.849 

.2321 

8.882 

.4795 

1.751 

71.27 

The  first  parameter  study  made  with  the  code  was  the  effect  of  pressure  on  the 
electron  number  density  as  shown  in  Pig.  2.  The  lower  curve  for  Cn  =  0  repre¬ 
sents  the  conditions  studied  in  Section  B  and  the  maximum  in  the  curve  that 
was  predicted  by  Professor  Rose  is  very  close,  indeed. 

Also  included  are  similar  curves  for  =  0.5  and  C,,  =1.0.  The  curve  for 
Cil  =  0.5  represents  the  additional  source  of  neon  metastable  states  equal  to 
1/2  the  neon  molecular  ion  recombination  rate.  The  curve  for  C-]_2_  =  1.0  repre¬ 
sents  the  limiting  case  where  all  of  the  neon  molecular  ion  recombination  rate 
goes  to  produce  neon  metastable  states.  As  shown  by  the  curves,  the  presence 
of  this  source  of  neon  metastable  states  would  greatly  enhance  the  electron 
number  density.  Very  little  is  known  about  this  reaction  so  that  a  reasonable 
value  of  is  very  difficult  to  come  by.  However,  with  our  new  understanding 
of  the  ion  generation  rate  in  neon  and  neonrargon  mixtures  and  the  forthcoming 
inpile  microwave  experiments  to  measure  the  electron  number  density  directly, 
it  may  be  possible  to  infer  a  value  of  this  quantity  since  the  rest  of  the 
important  reaction  kinetic  constants  are  reasonably  well  known . 

The  code  solution  also  yields  values  of  all  the  other  dependent  variables 

^2+,  Njj^^  Ny,  A+  and  A24.  and  values  of  these  variables  are  summarized  in 
Table  3  for  selected  runs . 

Effect  of  the  Concentration  of  the  Minor  Species,  Argon 

The  second  maximum  predicted  in  Section  B  was  that  the  electron  number 
density  would  peak  at ’the  ratio  of  argon/ neon  of  Aq/N^  =  1  x  10“^.  An  exhaus¬ 
tive  variation  of  this  variable  at  various  pressures  (N  )  was  not  made  but 
three  series  of  runs  of  Ng  versus  N^  were  made  for  Ao/No=  1  x  10" 3,  Aq/Nq  = 

1  X  10“  ,  and  A^/N^  =  1  x  10"5.  The  curves  are  presented  in  Fig.  3  for  the 
case  of  Oil  =  0.5  and  the  Qurve  for  Aq/Nq- 10“^  does  indeed  peak  at  the  high¬ 
est  electron  number  density.  The  peak  for  the  Aq/n^  _=  10’3  curve  has  been 
shifted  to  50  Torr  while  the  peak  for  the  A^/Nq  =  10-5  curve  has  been  shifted 
to  about  214-0  Torr . 
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Figure  4 


Effect  of  Electrode  Spacing 


As  pointed  out  previously  the  decrease  in  electron  number  density  with  decreas¬ 
ing  pressure  (at  low  pressures)  is  due  to  the  increased  diffusion  losses  of 
both  the  ions  and  metastable  states  (D^  =  K./No) .  .Now  the  diffusion  loss  (D^N^/ 
)  rate  also  depends  strongly  on  the  electrode  Spacing  (A  =  d/n)  so  that  it 
was  decided  to  also  make  a  preliminary  study  of  this  variable  particularly  since 
the  electron  transport  problem  for  an  energy  converter  implies  minimizing  the 
plasma  resistance  and  therefore  the  spacing . 

As  a  first  test  hhe  spacing  was  reduced  to  what  is  considered  a  lower  bound  from 
construction  considerations  of  .OSO  in.  or  A  =  .0162.  The  curve  of  Ng  vs  Nq 

reduced  A  =  0*5  is  shown  in  Fig.  4  together  with  the  curve  for 

A  =  0.1  for  comparison.  As  can  be  seen  Ng  is  greatly  reduced  and  the  maximum 
is  shifted  from  110  Torr  to  a  value  greater  than  400  Torr.  These  data  confirm 
that  the  optimization  on  pressure,  Aq/Nq,  and  geometry  must  also  include  the 
effects  on  plasma  resistivity  and  converter  performance . 
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SECTION  D 


PRODUCTION,  DimJSION  AND  CONVERSION  RATES  OF  IONS, 

METASTABLE  AND  EXCITED  STATES  IN  PENNING-TYPE  NEON-ARGON  PLASMAS 

ABSTRACT 

A  detailed  study  of  the  reaction  kinetics  of  a  fission  fragment  generated 
plasma  requires  knowledge  of  about  3O  reaction  rates  if  the  plasma  is 
created  using  a  Penning -type  gas  mixture .  These  reaction  rates  represent 
such  processes  as  diffusion,  recombination,  molecular-ion  formation,  charge 
exchange,  volume  ionization  and  destruction  in  the  gas  of  various  ionic, 
excited  and  metastable  states.  It  is  shown  that,  with  few  exceptions, 
these  reaction  rates  are  directly  known,  may  be  readily  extracted  from 
the  literature  or  may  be  estimated,  for  the  neon-argon  mix¬ 
ture  .  The  notable  exception  for  which  it  has  not  been  found  possible  to 
give  an  estimate,  is  the  production  rate  of  metastable  states  from  the  dis¬ 
sociative  recombination  of  a  molecular  ion  and  electron.  After  discussion 
of  all  the  reaction  rates,  the  data  are  collected  at  the  end  in  tabular 
form  suitable  for  reference  . 
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OBJECT 


The  purpose  of  this  work  was  to  make  a  compilation  of  the  fundamental  reac- 
tioa  rates  between  ions^  excited  states  and  neutral  atoms  in  a  r.eon-argon 
plasma  . 


SUMMARY 


Despite  the  large  number  of  reactions  which  can  occur  simultaneously  in  a 
neon/argon  plasma ^  sufficient  basic  data  have  been  compiled  for  worthwhile 
computational  studies  to  be  made  of  the  reaction  kinetics  of  this  system  . 


I .  INTRODUCTION 


Dense  plasmas  are  of  interest  in  the  field  of  thermionic  energy  conversion 
because  such  plasmas  can  be  used  to  neutralize  the  retarding  space  charge 
of  the  thermionic  electrons .  The  number  density  in  any  active  plasma  is  a 
function  of  both  the  ion  generation  rate  and  the  various  ion  loss  rates,  so 
that  in  order  to  achieve  maximum  number  density  it  is  necessary  to  study  in 
detail  the  dependence  of  the  ion  generation  and  loss  rates  on  such  parameters 
as  gas  species,  gas  density  and  electrode  configuration.  Very  recently  such 
a  study  of  reaction  kinetics  in  a  fission  fragment  generated  plasma  was  initi¬ 
ated  by  Professor  D.  J.  Rose  (l)  and  particle  density  equations  were  formu¬ 
lated  for  a  Penning-type  gas  mixture  which  are  applicable  to  a  neon-argon 
plasma  of  density  Ng  lO^^  electrons  cm“3 .  The  purpose  of  the  present  work 
was  to  make  a  compilation  of  the  various  reaction  rate  coefficients  in  neon- 
argon  which  were  needed  as  input  data  for  these  e'-uations  .  Table  I  summa¬ 
rizes  the  reactions  proposed  by  Professor  Rose  avid  will  be  continually  refer¬ 
red  to  in  this  report . 


n  .  GENERATION  RATE  OF  ATCMEC  IONS 
AND  EXCITED  STATES  BY  FISSION  FRAGMENTS 


Since  the  concentration  of  neutral  neon  atoms  N^  in  the  Penning  mixture  is 
much  larger  than  that  of  the  argon  atoms  Aq,  the  fission  fragments  in  passing 
through  the  gas  mainly  generate  atomic  neon  ions  N+,  metastable  excited  neon 
states  N^,  and  excited  neon  states  N^^,  that  are  not  metastable.  As  a  result 
of  excitation  transfer  and  the  imprisonment  of  resonance  radiation,  some  of 
the  Nx  states  decay  into  N^j^  states,  but  the  net  effect  of  such  transitions 
in  the  steady-state  will  be  to  increase  the  value  of  the  metastable  genera¬ 
tion  rate  S2 .  The  remairing  N^  states  will  either  (i)  undergo  radiative 
transitions  emitting  photons  which  escape  from  the  plasma  or  (ii)  produce 
N2+  molecular  ions  in  two  body  collisions  with  N^  (Hornbeck-Molnar  process ^2)] 
provided  the  Nx  state  has  an  energy  lying  within  the  range (3)  —  20.6  ev  to 
{=  21.6  ev)  where  E^^  is  ^he  first  ionization  energy  for  neon.  The 
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contribution  of  process (ii)  to  the  total  ion  generation  rate,  however,  cannot 
be  large  for  the  following  reasons.  Studies  by  Platzmanv^-'')  and  Utterback  and 
Millerw)  on  CC-particie  and  fission-fragment  ionization  of  noble  gases  indi¬ 
cate  that  the  value  of  +  S2)/Sg  is  i  .e  .  the  total  production  rate 

of  all  excited  states  is  about  one-half  the  ion  production  rate.  Furthermore, 
much  of  the  excitation  energy  will  be  distributed  among  excited  states  having 
energy  ^  20  .6  ev  and  which  are  thus  incapable  of  yielding  molecular  ions  via 
the  Hornbeck-Molnar  process.  The  production  rate  of  metastable  states  S2 
arising  from  fission  fragment  excitation  is  therefore  taken  to  be  much  larger 
than  the  production  rate  of  highly  excited  states  which  lead  to  molecular 
ion  formation,  i  .e  .  S2  -  S2_  ^  • 

The  value  of  may  now  be  determined  from  values  of  S  computed  for  pure  neon 
and  for  a  neon-argon  mixtui’e  (O.l7o  Ar)  under  the  same  conditions.  At  a  pres¬ 
sure  of  2k0  Torr,  SQ(Ng)^o  =  2.22  x  lOl^  cm“3  sec"l  immediately  a.djacent  to 
the  uranium  surface  for  a  neutron  flux  value  of  0  =  O.99  x  10l3  cm-2  sec"l. 

For  the  same  values  of  pressure  and  neutron  flux,  the  total  ion  generation 
rate  for  the  neon-argon  mixture  is  Sg(NeAr) .N^  =  3 .13  x  10^^  cm'J  sec"^.  Neg¬ 
lecting  metastable -metastable  collisions,  (the  density  of  metastable  states 
is  negligible  compared  with  the  density  of  neutral  atoms),  and  considering 
the  sheath  above  the  uranium,  where  recombination  of  ions  does  not  occur  and 
cannot  thus  contribute  to  metastable  production,  the  conservation  equation 
for  metastable  states  is 

\  \  No  N^  .  Ao  qi  J 

D  N 

^  ^eff  \  ^m  "^i 

2  A  2 

Replacing  y  t)y  -  l/iY  ,  the  reciprocal  squared  of  the  diffusion  length,  is 
valid  only  if  diffusion  loss  dominates.  However  if  the  diffusion  loss  is  small 
compared  with  the  other  loss  terms,  then  its  exact  value  is  of  no  consequence 
in  the  above  equation,  (3.6  x  10-15  cm3  sec"^)  is  an  effective  2-body 

destruction  frequency  for  neon  metastable  states,  q^^  (2.5  x  lO-l^  cm^)  is  the 
cross-section  for  ionization  of  argon  by  a  neon  metastable  and  v  (6 .9  x  10^ 
cm  sec“l)  is  the  average  velocity  of  approach  of  these  two  particles;  D  = 

0.65  cm2  sec“l  for  Kq  =  8.5  x  IQIS  cm'5,A^  =  9  x  10-3  cm^  and  A^  =  8.5^  10l5 
cm"3.  All  these  nuir.erical  values  are  discussed  in  detail  below.  Substituting 
the  niunerical  values  in  Eq.(l)  yields 

^2^0  ^  72  N^  +  3.1  X  lo'^  N^  +  1.5  X  10^ 

=  1.8  X  10^  Njji.  (2) 

The  ion  generation  rate  in  neon-argon  consists  of  the  ion  generation  rate  in 
pure  neon  plus  the  generation  rate  of  argon  ions  resulting  from  collisions  of 
neutral  argon  atoms  with  neon  metastable  states,  i  .e . 
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^_(NeAr)*  N  =  S_(Ne).  N  +  N  A  q.  v 
U  o  O'  o  m  o 

(3.13  -  2.22)  X  10^^  =  1.5  X  10^  N 

m 

'  N  =  6 .1  X  10^*^  cm"^  (3) 

Using  Eqs  .  (2)  and  (3)^  the  generation  rate  of  neon  metastable  states  is  found 
to  be 

=  (1.8  X  10^)(6.1  X  10^^) 

^  ^ - 

8.5  X  10^° 

-  1.3  X  10  ^  sec~^ 

III  .  DIFFUSION  COEFFICIENTS  FOR 
KETASTABLES,  ATOMIC  AUD  MOLECULAR  IONS 


The ^arabipolar  diffusion  coefficients  for  Ne  ,  NOp"^,  Ar"*"  and  At 2^  in  neon  at 
300°K  are  determined  from  measurements  of  ion  mobility  Pp’^Cat  standard  condi¬ 
tion  of  273 °K  and  76O  Torr)  using  the  following  relationships: 

+ 

- =  .026  ev  for  T  =  300°K,  and 

M-q  ^ 

=  2D  for  thermal  equilibrium  =  T" . 

The  ambipolar  diffusion  coefficient  at  unit  atom  density,  is  defined  by 
^a“^o'^a  ^q(2  .69  X  10  ^  cm~3)  ±s  the  number  density  of  atoms  at  760  Torr 

and  273 "K.  Therefore  K  is  given  by 

3. 

=  2  X  .026  X  2.69  X  10^^  =  1.4  X  10^®  cm"^  sec'^ 

2  1  1 

(for  q^  in  cm  V  sec'  )  (4) 

K  (Ne"^)  in  neon 

J-;?- _ 

f  2  1 

Uo(Ne  )  =  4  .1  cm  (v*sec)  from  the. work  of  Oskam  and.  Mit- 

telstadt'^^  and  McDaniel'^'^ . 

Therefore  using  Eq .  (4),  is  given  by 

K,  =  5 *7  X  10^^  cm  ^  sec'^ 

A.  y3. 
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[io(ife2^)  =6.5  cm^(v.sec)"^  frow  Oskam  and  Mittelstadt^^^ 

and  McDaniel(7) 

'  Kq  ^  =  9*1  X  10^^  cm  sec”^ 

K  (Ar^)  in  neon 

j}^ _ _ _ 

The  value  of  ho(^yo)nF  =  7-^  cm^(v.sec)“^  is  taken  from  the  data  of 

Chanin^and  BiondlVO;9J  showing  mobility  as  a  function  of  ion  mass  in  neon. 
Primarily  Chanin  and  Biondi  determined  the  mobility  of  Hg”^  ions  in  neon,  but 
the  results  are  in  agreement  with  the  extrapolation  of  mobility  versus  mass 
niMber  curves  obtained  for  the  alkali  ions  in  the  noble  gas  by  Tyndall , 

These  data  are  considered  appropriate  for  p,  (Ar'*')  in  neon  since  resonant 
charge  transfer  effects  (important  for  an  atomic  ion  diffusing  in  its  parent 
atomic  gas)  are  absent. 


'  K  =  1 .0  X  10^^  cm  ^  sec'^ 

«  • 

^4  a 

Po(Ar2^)  in  neon  =6.8  cm^(v.sec)’^  from  Chanin  and 

Biondi^ 9) 

^4  a  ”  9*5  X  lO"*"  ■  cm"^  sec~^ 

The  ^diffusion  coefficient  for  the  neon  metastable  states  3P2  and  3Pq  at 
300  K  has  been  determined  by  Phelps from  optical  absorption  measurements 
on  excited  neon  atoms  following  a  pulsed  discharge.  For  both  states,  it  was 
found  that 

V  c;  c:  v  1  “7 

=  5  *5  X  10  cm  sec 

The  characteristic  diffusion  length.-ZVfor  infinite  parallel  plate  geometry 
and  for  the  lowest  diffusion  mode'9/  is  (L/n)  where  L  is  the  plate  separation. 
For  L  =  0.3  cm, 

=  9  X  10”^  cm^ 


IV.  ELECTRON-ION  RECOMBINATION  COEFITCIENTS 
FOR  NEON  AiJD  ARGON  ATOMIC  AilD  MOLECULAR  IONS 


Tliere  is  a  dearth  of  reliable  experimental  determinations  of  atomic  ion- 
electron  recombination  coefficients(Ci,Cp,CT,C-,n  C-in)  „ 

"thG  fact,  that  except  for  dense  plasrnas  10^  ern"-^)  these  processes  are 

normally  small  in  comparison  vrLth  other  loss  processes  and  are  consequently 
difficult  to  measure.  The  process  of  radiative  recombination(92_,Ci3T 
ever,  been  Investigated  theoretically  in  some  detail  for  hydrogenic  ions 
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(Bates  and  Dalgarno^  and  "bhesc  data  are  considered  to  be  of  high  accuracy. 
For  the  atomic  hydrogen  ion  the  total  radiative  recombination  coefficient 
at  250°K  is  4.84  x  10“^cm^  sec"^  and  at  5(^0°K  is  3.12  x  10“i2  c^3  sec"^. 
Approximate  values  of  the_^total  radiative  recombination  coefficient  for  more 
compD.ex  ions  (He'*’,  Li*",  C  ,  N^,  O'*",  Ne'*',  ,  KT*’ )  have  been  obtained  at  250°K 

on  the  assumption  that  all  excited  levels  are  hydrogenic  and  using  knovm  val¬ 
ues  of  the  threshold  photoionization  cross  sections  from  the  ground  levels 
of  the  atoms  concerned.  For  Ne*"  at  250°K,  the  radiative  recombination  is 


and  for  the  remainder  of  the  Ions  listed  above,  the 


3  .4  X  10"^^  cm^  sec“^ 

recombination  coefficient  ranges  from  3.0  x  10-12  to  4 .8  x  lO"!^  cm'  sec 
showing  that  the  values  do  not  differ  much  from  one  singly  charged  ion  to 
another.  For  He'*'  at  300°K,  it  is  considered  that  the  value  of  is  close 
to  3  .2  X  10“12  cml  sec"l  and  in  the  absence  of  explicit  data  for  Ar”^,  Cp3  is 
taken  4  x  10  cm^  sec-1  by  analogy  with  the  atomic  ions  considered  above. 


3  *2  X  10  cm^  sec~^ 

C^g  ^  4  X  10  cm^  sec"^ 


Cp,C^  and  C.g  represent  3-‘body  atomic  ion-electron  recombination  coefficients, 
where  the  third  body  for  Cp  and  is  a  neutral  neon  atom,  and  for  is  an 
electron.  In  this  radiationless  recombination,  the  function  of  the  third  body 
is  to  remove  the  surplus  energy  released  by  recombination.  For  plasma  densi¬ 
ties  Ng  —  10l2  cm"3  and  gas  pressurp  >  30  Torr,  Ng/N©  ^  10"°  so  that  if  C-jNg 
is  to  be *^2^0^  ^3/^2  ^^^st  be  ^10^.  It  is  inconceivable  that  an  electron 
would  be  a  million  times  more  efficient  than  a  gas  atom  for  the  removal  of  sur¬ 


plus  energy,  so  that  under  the  present  conditions  C-^<<  C2  N^/Ng  and  need  not  be 
considered  further  .  (‘Ihe  coefficient  corresponding  to  for  A'*'  would  also 
lead  to  a  negligible  recombination  loss  in  comparison  with  C;]^^Ng  and  has  there¬ 
fore  been  omitted  from  the  equation  for  A"*”,  Table  I.) 


Approximate  values  of  Cp  and  may  be  obtained  from  an  extension  of  the 
classical  theory  of  Thomson  (Massey  and  Burhop^^S)).  The  general  expression 
for  C2  and  is 


(8kT/n  m)^^^  8  n  r  m 

0=— - 

3  M  q^ 


(5) 


where  m  is  the  electron  mass,  M  is  the  mass  of  the  gas  atom,  n  is  the  number 
density  of  atoms  at  1  Torr  and  273°K  and  qg  is  the  total  elastic  cross  section 
for  electron-atom  collisions  .  He combi nation  is  assumed  to  occur  when  electrons 
and  ions  approach  each  other  within  a  distance  r^  defined  by 


r 

o 


(C) 


=  3-71  X  10 


-6 


cm  at  T  =  300 °K. 
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On  this  simplified  viev/  oJ  3-tody  recombination,  the  important  consideration 
i.j  the  rate  at  which  electron::  can  lose  energy  in  collisions  with  neutral 
atoms,  v:hlch  is  dependent  on  M  and  for  a  given  temperature.  Since  >  >  A  , 
M  and  are  respectively  the  mass  of  the  neon  atom  and  the  neon  atom-electron 
cross  section  in  the  calculation  for  C2  and  also  for  Cqq  so  that  Cp  =  Cqn .  The 
value  of  for  neon  is  deteraiTied  from  the  experimental  values  of’P  ,  the 
probability  of  collision  pF  electrons  in  noon  (q  =  0.283  x  10^^  P^)^obtained 

by  Gilardini  and  Brovm.^1^^  At  300°K,  P^  =  1 .  ,(c:;i  .Torr)-l,  so  that  q  = 

3.7  X  10-1^  cm2.  c  ae 

"  (8  X  1.38  X  10~^'^  X  300/3.14  X  0.11  X  10~^'^)  ^ 

3  X  (20  X  1.87  X  10  "■')  X  3  .7  X  l  X  ('^,5)4  X  1'^^''^)  ' 

X  8  X  3-1't  X  (3.71  X  10”^)  X  ^.11  X  10"^^ 


Cp;  ~  2 .7  X  10  cm'^  sec  ^ 

In  contrast  to  the  dearth  ol‘  experimental  determinations  cf  atomic  ion- 
electron  recombination  coefficients,  many  reliable  experimental  values  are 
available  for  molecular  ion-electron  dissociative  recombination  (Cqg,  Cp-|  )  . 
This  process  can  occui*  ra]vidly  and  is  often  the  dominant  loss  mechanism  in 
plasmas  for  which  25  5  x  lO^T  cm'^  and  10l2  cm'^’ .  The  latest  pub- 

lished  values  ^  of  and  Cpq  at  300°K  given  below  are  in  good  agreement 
with  earlier  values  obtained  by  BiondiUo,'  and  Biondi  and  Brow.^^^^ 


-  2.2  X  10  cm’  sec  ^  from  Oskam  and  Mittelstadt ^ 
*"21  "  X  10  ^  cm'*  sec  ^  from  Oskam  and  Mittelstadt ^ 


V .  PRODUCTION  OF 

EXCITED  STATES  FROM  RECainEMTEON  PROCESGES 


elec- 

O 


imi- 


It  is  energetically  possible  that  the  recombination  of  neon  ions  and 
trons  can  lead  to  the  production  of  excited  neon  states  (C3,C3_q,C-| -.  ) 
lai  coefficients  for  the  production  of  excited  ai'gon  states  are  not  con¬ 
sidered,  since  such  states,  cannot  usefully  contribute  by  any  scheme  to  the 
total  ion  density.  Tlieoretical  studies  by  Bates  and  Daigarno(12)  radia¬ 
tive  recomljination  have  ;diov.Ti  that  in  the  case  of  ions,  electron  capture 
to  tlie  ground  s.tate  i::  12  times  more  probable  than  electron  capture  to  a 

liighly  exciood  state.  On  this  basis  Cp  for  radiative  recombination  (Cq)  is 

in  fact,  be  a  different  coefficient  in  Table  1  for  the 
'tator',  from  the  Tfiomson  3-Eody  recombination  process  (Cp 
for  efficient  removal  of  the  energy  of 
of  tlie  atoms  being  in  an  excited  state 
small.  Fuj-thenaore  since  there  is  no 


0.1„  Tlma-e  should, 
production  of’  excited  ; 

However,  since  3  bodier:  arc  available 
recombination,  the  probability  f  one 
after  recombination  is  expected  to  be 


) 


published  work  on  the  production  o:f  excited  states  by  3-body  recombination. 
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and  more  importantly,  since  atomic  ion-electron  recombination  is  not  a  sig- 
ni  leant  loss  process  for  neon  under  the  present  conditions,  the  separation 
01  C3  into  two  coefficients  is  not  necessary.  Similar  arguments  apply  to 
the  production  of  metastable  states  C]_q  from  atomic  ion-electron  recombina¬ 
tion.  Since  the  energy  of  recombination  has  to  be  shared  among  the  and 
^  states,  it  is  expected  that 0.1  would  be  an  upper  bound  for  both  Cp  and 
Cio.  but  much  significance  should  not  be  attached  to  this  value . 

"^8^  "10  < 

Contrary  to  the  unimportance  of  Cq  and  C-j_Q,  precise  knowledge  of  C-i-i  is  highly 
relevant  to  the  present  study  because  of  the  large  value  of  C-,g  and  the  con¬ 
sequent  possibility  of  a  significant  source  of  metastable  states  via  C  C  ^ 
^2+^e  ■  Unfortunately,  little  is  known  about  the  production  of  excited  or^ 
metastable  states  from  dissociative  recombination,  although  it  is  generally 
accepted  that  the  dissociative  recombination  of  diatomic  noble  gas  ions  leads 
to  dissociative  atoms  being  left  in  an  excited  state  (Biondi  and 

Holstein  ).  Recent  work  by  Biondi  et  al .  on  the  dissociative  recombina¬ 
tion  mechanism  in  neon'^  and  helium(20)  plasmas  using  microwave,  optical 
spectrometric  and  interferometric  techniques,  has  indicated  that  dissocia¬ 
tive  recombination  into  high-lying  atomic  excited  states  is  a  likely  process . 
Since  the  transition  probabilities  of  some  of  the  exc ited  states  (above  the 
metastable  states)  to  the  ground  state  is  finite,  C,-.  is  likely  to  be  less 
than  one,  but  further  identification  of  the  high-lying  excited  states  is 
required  for  estimations  of  . 


qi<  1 

VI .  THREE -BODY  MOLECbTLAR  ION  FORMATION 


The  coefficient  for  the  conversion  of  atomic  neon  ions  into  molecular  ions 
in  3-body  collisions  involving  the  neon  ion  and  two  neon  atoms,  has  been 
determined  by  Beaty  and  Patterson.^^^^ 

Cl  =  5 *8  X  10  cm^  sec"^ 


The  coefficient  C^,  representing  a  collision  between  an  atomic  neon  ion,  a 
neon  atom  and  an  argon  atom,  is  not  known,  but  since  the  argon  takes  no  part 
in  the  collision  except  to  satisfy  conditions  for  conservation  of  momentum 
will  be  approximately  equal  to  .  More  precise  knowledge  of^Cc  is  not^ 
required  for  the  present  conditions  since  Nq>>  and  thus  C^Nq"^  >>  C^N^A^ . 

5  X  10"^^  cm^  sec"^ 

^e  production  rate  of  molecular  argon  ions,  C20,  from  3-body  collisions 
involving  an  atomic  argon  ion,  an  argon  atom  and  a  neon  atom  has  been  inves¬ 
tigated  experimentally  by  0skam.v22;  The  relevant  conversion  frequenev  is 
found  to  be  j  • 

=  2.5  X  10'3^  cm^  sec'^ 


8 


vn  .  WO-BODY  CJIARGE  EXCHANGE 


Two-body  charge  exchange  (C5,C2_y)  ir  energetically  possible  when  the  poten¬ 
tial  energy  of  the  incident  ion  is  greater  than  or  equal  to  the  ionization 
energy  of  the  target  atom.  Massey  and  Burhop^^^)  pave  indicated  that  the 
cross  section  for  charge  exchange  will  be  maximum  when  the  collision  time 
a/v  is  approximately  equal  to  the  time  for  electronic  transition  h//.E.  The 
quantity  a  is  of  order  10~8  cm,  the  diameter  of  an  atom,  v  is  the  relative 
velocity  of  approach  of  the  two  particles,  h  is  Planck's  constant  and  Ae  is 
the  difference  in  the  ionization  energies  of  the  unlike  colliding  particles . 
For  maximum  cross  section 


V 

m 


aAE 

h 


(7) 


For  the  reaction  represented  by  C^(Ne^ Ar”^) ,  Ae,-.6  ev;  for  (Ne^^ — >Ar^), 

5  ev  from  the  data  of  Mason  and  Vanderslice . ^ ^ ^  Using  these  values,  v 
becomes  respectively  =  10 '  cm  sec“l  and  ^  8  x  10^  cm  sec"!.  For  temperature  T°K, 

V  =  (8kT(Mq  +  M^)/n  where  M^^  and  M2_represent  the  mass  of  particles  1 

and  2  involved  in  the  collision.  At  300°K,  v  is  given  by 


v(Ne,  Ar)  = 


8  X  1.38  X  10"^^  X  300  X  60 


3 .14  X  20  X  4o  X  1 .66  X  10 


4  -1 

=  6 .9  X  10  cm  sec 


Similarly  v  (^62^  Ar) 


4  -1 

5 .7  X  10  cm  sec 


Therefore,  since  v  «  v^  for  both  the  atomic  and  molecular  neon  ion  it  is  evi¬ 
dent  that  the  charge  exchange  cross  sections  fq^x)  will  be  small  (compared  to 
gas  kinetic  q  10  cm  ^  for  the  values  of  A  E  and  low  relative  velocities 
considered  here,  q^^  (Ne  ^Ar^)  at  300°K  can  be  approximately  estimated  by 
extrapolation  of  the  theoretical  data  obtained  by  Hasted for  energies 
~  12  to  30  Gv,  but  no  data  are  available  for  the  molecular  neon  ion.  How¬ 
ever  ,  despite  the  difference  in  value  for  A  E  for  the  atomic  and  molecular 
ion,  the  charge  exchange  cross  sections  at  300°K  (far  from  ^he  maximum  cross 
section)  will  tend  to  converge  towards  zero  so  that  q^^  (Ne--->Ar+)  is  taken 
to  be  q„,,(Ne — ^Ar"*").  From  Har,tpdv23).  extrapolated  value  for  the  prob- 

Therefore,  since 


±QX'^  /  •  From  Hasted , 

ability  o-r  charge  exchange  (Ne”*" - >Ar'^)  cm"^ 


,(Ne^^ — > Ar"*")  . 

aharge  exchan.  , 

^6^^17  given  by  v  O.283  x  IQ-l^^ 


^6 ’^17 


3.9  X  10'^^  cm^  sec'^ 


VIII .  COLLISIONS  INVOLVING  METASTABLE  STATES 


The  concentration  Nn;  of  neon  metastable  states  in  the  present  study  is  suf¬ 
ficiently  low  (^  lol^-lO^l  cm"^.  Section  II )  that  collisions  between 
pairs  of  metastable  states  will  be  very  small  in  comparison  with  collisions 
between  a  metastable  state  and  a  neutral  atom.  Even  so,  an  order-of -magnitude 
estimate  of  the  coefficients  C.^,C^2  useful  in  establishing  that  the 


9 
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of  mctastiable -metias'table  collisions  is  negligibly  small .  The 
metastable-metastable  ionizing  cross  section  does  not  appear  to  have  been 
determined  for  neon,  but  the  cross  section  has  been  estimated  by  Phelps'^^^ 
for  the  case  of  helium.  Phelps  obtains  10"^^  cm^,  and  on  the  basis  that 

the  ionizing  cross  section  for  neon  metastables  is  of  similar  order  of  magni¬ 
tude,  Cj(=  v(Ne,Ne)  where  v  =  8  x  10  cm  sec"^)  becomes 

^  8  X  10  cm^  sec  ^ 

( 25 ) 

In  some  recent  work,  Myers^  '  has  suggested  that  metastable-metastable 
ionizing  collisions  lead  to  the  production  of  molecular  ions  and  not  atomic 
ions.  However,  since  studies  presently  being  conducted  by  Biondi,(20) 
together  with  earlier  work} 26, 27)  strongly  support  the  conclusion  that  meta¬ 
stable  -metastable  ionizing  collisions  yield  atomic  ions,  the  possible  pro¬ 
duction  of  molecular  ions  from  metastables  is  not  considered  here. 

All  of  the  collisions  between  metastable  states  do  not  lead  to  atomic  ion 
production  so  that  C-j^^  ^  although  not  necessarily  by  a  significant  amount . 
Because  of  the  large  size"  of  the  metastable  atom,  the  cross  section  for 
collision  between  metastable  states  will  be  larger  than  the  gas  kinetic  cross 
section  for  collision  of  neutral  states.  However,  the  value  of  q^  above  is 
already  so  large  that  it  suggests  that  almost  all  of  the  metastable-metastable 
collisions  lead  to  iorl’ation.  For  the  present  purposes  such  an  estimate  is 
considered  sufficieni  . 

C12  8  X  10”^^  cm^  sec”^ 

The  predominant  destruction  processes  for  neon  metastables  in  the  volume  of 
the  gas  are  2-body  collisions  (C^.)  with  neon  atoms,  3-body  collisions  (C3_4) 
with  neon  atoms,  and  2-body  ionizing  collisions  with  argon  atoms  (Pen¬ 

ning  effect).  All  these  processes  have  been  experimentally  Investigated  in 
some  detail.  The  cross  section,  q^  for  ionization  of  an  argon  atom  by  a  neon 
metastable,  has  been  determined  by  Biondi^^°^  as  2.6  x  10-16  cm^ .  The  product 
of  and  v(lfe,A)  yields  . 

=  1.8  X  10  cm^  sec'^. 

Low  pressure  microwave  studies  (p  <  2.5  Torr)  by  Biondi^^®^  of  the  destruc¬ 
tion  of  neon  metastables  in  neon  indicated  that  the  predominant  volume  loss 
was  by  2-body  collisions  with  a  destruction  cross  section  q^  =  8.9  x  10"^6 
cm2.  On  this  basis  =  7  •!  x  10"^5  cm3  However,  subsequent  more 

detailed  studies  were  undertaken  by  Phelps''"^"^  ^  using  optical  absorption  tech¬ 
niques  to  measure  decay  rates  of  excited  atoms,  and  these  investigations 
showed  that  at  higher  neon  densities  (lolT_iol9  atoms  cm"3)  the  volume 
destruction  of  metastable  states  was  complex,  involving  excitation  and  de-excl- 
tation  of  the  3P^  and  3P2  metastable  states  to  the  nearby  3Pj_  resonant  radi¬ 
ating  state  and  also  involving  3-body  destruction.  The  3-body  destruction 
term  became  increasingly  important  as  expected  for  Nq  >  10l°  atom  cm-3,  and 
for  No  =  I0I9  atom  cm"3  the  destruction  rate  was  given  by  Cpk  where 
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Ci4  =  5  X  10”^ cm  nec”  at  300°K.  Experimental  valuer,  of  the  decay  rate 
of  metastabler  were  obtained  for  5  ^  10^5  IqI*^'*  atoms  cm"  "i  and  Phelps 

was  successfully  able  to  fit  all  the  experimental  points  with  a  single  expies- 
sion  which  contained  diffusion  loss  dominant  at  low  pressures)  1-body 

destruction  loss  Nq  dominant  at  high  pressures)  and  a  tenn  which  took 

into  account  the  excitation  and  de-excitation  transfers  which  occur  between  the 
3P  ,  iPg  and  3P^  states  .  Use  of  the  latter  teriTi  requires  a  knowledge  of  the 
radios  of  the  number  densities  of  both  the  3^0  states  to 

the  3P2^  state  and  these  data  have  been  given  by  Phelps  for  the  range  of  values 
of  Nq  investigated  . 

For  the  present  studies,  however,  Phelps' full  expression  may  be  considerably 
simplified  in  the  following  manner,  without  significant  loss  of  acc\rracy . 

The  characteristic  diffusion  length  used  by  Phelps  was  given  by  cm^ 

so  that  in  comparison  with  the  present  case  (A^  "  *^*0^9  cm^.  Section  III); 
diffusion  loss  was  much  smaller.  For  values  of  A^<.0.01  cm  ,  diffusion  loss 
will  dominate  for  N^,;S  4  x  10^'  atoniQcm‘A  For  the  range  4  x  IOIT.S 
8.5  X  10^  atom  cm"'>  (Nq  =  8.5  x  10^°  is  p  =  240  Torr)  the  total  volume 
destruction  frequency  given  by  Phelps  may  be  approximated  by  an  effective 
2-body  destruction  cross  section  =04*5  x  10"^^  cm^  or  Vgpp  -  3 -6  x  10-  5 

cm^  sec“^  where  VgppN^  =  i-  Nq^.  A  measure  of  the  degree  of  approxi¬ 

mation  is  given  in  Table  H  . 


TABLE  II 

Comparison  of  Effective  Decay  Frequency  for 
2-Body  Metastable  Destruction  with  Phelps'  Experimental  Values 


N  (atom 
0 

cm”^) 

Effective  Decay  Frequency  for 
2-body  Destructive  Collision, 
VeffNoSec"^  with  v^^^=3'6xlO" 
cm8  sec" 

15 

Phelps'  Experi¬ 
mental  Value  of 
Total  Decay  Fre¬ 
quency 

8  .5  X 

10^® 

3-1 

-1 

X  10  sec 

4.0 

X 

■,  -1 

10  sec 

4  X 

1018 

1 .4 

X  10^ 

1.2 

X 

10^ 

1  X 

10^8 

3.6 

X  10^ 

2.8 

X 

10^ 

8  X 

10^’^ 

2.9 

X  lo"* 

2.5 

X 

10^ 

4  X 

10^’^ 

1.4 

X  lo"^ 

1.4 

X 

10^ 

Evidently, 

for  A^<  0.01 

cm^  and  4  . 

X  10^^ :5  ®' 

,5  X 

10^8 

atom 

cm"^,  C^^  + 

C14N0 

may  ' 

be  approximated  by  v^pp  where 

V  = 

eff 

3  .6  X  10"^^  cm'^ 

sec 

■1 

rX.  COLLISIONS  INVOLVING  HIGHLY- EXCITED  STATES 


It  has  been  mentioned  pi-eviously 
lar  neon  ions  from  highly-excited 


(Section  II )  that  the  production  of .molecu- 


U  X  w  1 1  xa.  /  Uiiciu  ^  A  vy  VA  04.V-.  1  V-/X 

states  Ny  (Hornbeck-Molnar  process''^^ 


IS 
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not  an  important  process  under  the  present  conditions .  However,  typical  val¬ 
ues  of  where  is  the  cross  section  for  ion  formation  from  Njr  and  t 

is  the  mean  lifetime  of  N^,  have  been  determined  by  Hornbeck^^^^  for  helium, 
neon  and  argon  and  the  value  for  neon  is  a#  5  x  10"^3  cm^  sec  .  This  value  of 
yields  the  ratio  of  the  number  of  molecular  ions  formed  sec-1  to  the 

number  of  excited  atoms  decaying  sec“l  under  given  conditions  since  this  ratio 
is 


N  N  q_  ^ 
X  o  ^ 

iTTr 

x' 


=  N  V 


(8) 


l8 

For  N  =  3.54  X  10  atoms  cm“3  (p  =  100  Torr),  v  a  t  =»  l4  which  indicates 
that  III  molecular  ions  will  be  formed  from  for  each  %  Ih^t  radiatively 
decays .  Evidently  for  the  situation  of  copious  production  of  excited  states 
with  sufficient  energy  for  the  production  of  Neo^  (within  1  ev  of 
Ei(Ng)  =  21.6  ev)^3J  this  process  would  be  a  significant  source  of  molecular 
ion  generation . 


Ib  is  worth  noting  that  if  q^  is  of  large  atomic  magnitude  ^10"^^  cm^  (since 
the  state  is  highly  excited),  then  t  5  x  10"8  gee  which  appears  reason¬ 
able.  The  coefficient  would  then  be  of  order  of  magnitude  8  x  10"^^  cm^sec"\ 


5  X  10’® 


sec 


8  X  10”^^  cm^  sec"^ 


X.  COLLEC'ECON  OF  DATA 


The  production,  diffusion  and  conversion  coefficients  of  ions  and  metastable 
states  discussed  in  this  report  are  summarized  for  convenience  in  Table  III. 
All  the  data  are  evaluated  for  a  temperature  of  300“K. 


12 


Process 


TABLE  ni 

Rate  Coefficients  for  Ne/Ar  Plasma  at  300°K 

Symbol 


Value 


Production  Ne  by  fission  fragments 

Diffusion  Ne^  in  Ne 

2-body  radiative  recombination  Ne^ 


o 


K 


l,a 


+ 


3 -body  Thomson  recombination  Ne 
3-body  2-electron  recombination  Ne 
3-body  Ne 
3 -body  Ne^ 

2-body  Ne 


Ne^"*^  with  2  Ne 


2 .6  X  10”^  sec~^ 

c;  -7  -1  -1 

5.7  X  10  cm  sec 

_  „  ,»-12  3  -1 

3.2  X  10  cm  sec 

^  „  ,.-30  6  -1 

2 .  f  X  10  cm  sec 

<<  C„N  /N  cm*^  sec”^ 

2  o  e  ^ 

5 .8  X  10”^  cm  sec” 


Ne, 

< 

->Ar‘ 


with  Ne,Ar 


m 


+ 


Metastable -metastable  Ne“ — >Ne 
Production  Ne  by  fission  fragments 
Production  Ne  by  atomic  recombination 
Mean  lifetime  Ne 

X 

2  body  ion  formation  Ne - >  Ne„ 

X  2 

Production  Ne^  by  fission  fragments 

Diffusion  Ne  in  Ne 
m 

Production  Ne^^  by  recombination  of  Ne 

Production  Ne  by  recombination  of  Neo' 
m  2 

Metastable -metastable  Ne’^  destruction  C 

Destruction  Ne  via  Ne 
m 

Penning  destruction  Ne  +Ar-^Ar 
+  ^ 

Diffusion  Ne2  in  Ne 


5 

X 

10”^^ 

6 

cm 

-1 

sec 

4 

X 

io”^3 

.n,3 

cm 

-1 

sec 

8 

X 

10-1° 

3 

cm 

-1 

sec 

-1 


'6 


T 

^2 

K 

m 

^10 
^11 
12 
^eff^ 


K 


'15 

.  2, a 

2-body  dissociative  recombination  Ne„ 

.  ^  2  16 


2-body  Ne^  — ^  Ar 
Diffusion  Ar”^  in  Ne 

2- body  radiative  recombination  Ar^ 

3- body  Thomson  recombination  Ar^ 
3-body  Ar"*"  +  Ar  +  Ne — 

Diffusion  Ar^  in  Ne 

2-body  dissociative  recombination  Ar^^ 
Characteristic  diffusion  length 


K 


'17 

3, a 
^18 
^19 
*^20 
^U,a 

% 

K 


«  Sg  sec 

0.1 

q  -I  p.-8 

5  X  10  sec 

8  X  10  cm'  sec 

1 .3  X  10”^  sec”^ 

5 -5  X  10^^  cm”^  sec"^ 

S.  0.1 

<1.0 

.  n  T^-IO  3  -1 

o  10  cm  sec 

=C  +C  I  N  )*  3*6  X  10  cm^  sec”^ 

1 .8  X  10  cm"^  sec 

9.1  X  10^^  cm”^  sec”^ 

2.2  X  10  cm^  sec'^ 

-13  3  -1 

cm  sec 

T  ^  ,a9  -1  -1 

1.0  X  10  cm  sec 

,  ,^-12  3  -1 

4  X  10  cm  sec 

2.7  X  10 

2.5  X  10 

9.5  X  10 


4  X  10 


9  X  10 


-3 


*The  approximation  of  _  +  C. ,  N  to  v  „„  is  valid  for  A  <  0.01  cm  and 
I*  X  lOlT  5.N  s8.5  X  MlS 


>-30  6 

-1 

>  cm 

sec 

,-31  6 

-1 

1  cm 

sec 

,18  -1 

-1 

>  cm 

sec 

rT  3 

-1 

sec 

1  2 

cm 

cm^  and 
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SECTION  E 


ION  GENE.RATION  RATE  MEASUREMENTS  IN  AN  ARGON-CESIUM  MIXTURE 


ABSTRACT 

Ion  generation  rate  measurements  were  made  on  an  argon-cesium  gas  mixture 
using  modified  versions  of  the  noble  gas  ionization  tube  .  Data  from  the  first 
two  tubes  were  inconclusive  due  to  the  difficulty  of  transporting  cesium  vapor 
in  argon  (at  2k0  torr)  for  the  tube  configuration  used.  A  study  of  cesium  trans¬ 
port  in  argon  showed  that  a  period  of  one  to  two  hours  is  needed  to  establish 
equilibriimi  and  this  suggested  further  modifications  to  tube  design.  The  third 
ion  tube,  thus  modified,  operated  successfully  in  the  reactor.  The  data  indi¬ 
cated  that  for^Cs/Ar  of  10  ^  at  an  argon  pressure  of  240  torr,  the  I-V  data 
followed  a  (v)  relationship  and  indicated  that  the  ion  generation  rate  is 
significantly  higher  than  that  for  Ne-Ar  or  the  pure  gases  for  the  same  condition. 
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1 .  HITRODUCTION 

The  results  for  Ne-Ar  (Gection  C)  suggested  thut  sorpe  other  mixed  gus  sys¬ 
tems  might  be  capable  of  sustaining  high  ion  densities  v^hen  ionized  by  fission 
fragments  .  •  For  applications  in  thermionic  converters  argon  seeded  vdth  ces.ium 
was  considered  to  be  an  attractive  Penning-type  sys.tem  because  of  the  small 
electron  scattering  cross  section  of  argon  and  the  lov^  ionisation  potential  of 
cesium . 

The  ion  generation  rate  program  was  thus  oriented  tov/ardc  a  study  of  Ar-Cs 
plasma  generated  by  fission  fragments.  For  this  purpose  the  noble  gas  ion  tube 
v;as  modified  to  be  compatible  vrith  cesium.  Since  currents  of  order  ma  .  or  less 
were  expected^  the  measuring  circuit  had  to  be  designed  to  nullify  the  influence 
of  leakage  current  across  the  ceramic  insulator . 

The  results  of  the  xirst  tv/o  tubes  are  described  briefly.  These  suggested 
that  the  results  were  influenced  by  atom  transport  effects  in  the  Ar-Cs  system. 
The  transport  of  cesium  in  argon  v/as  subsequently  studied  in  the  laboratory  and 
this  led  to  a  major  modification  in  the  ion  tube  geometry.  The  results  from 
the  data  obtained  inpile  vdth  this  tube  are  presented  in  terms  of  the  ion  trans¬ 
port  and  ion  generation  rate  models  developed  for  the  noble  gas  system. 

II .  CIRCUIT 

In  the  accurate  determination  of  small  (  ^  10  'amps)  ionization  currents  in 
diodes  filled  with  cesium  vapor,  care  must  be  taken  to  ensure  that  leakage  cur¬ 
rents  across  cesium-covered  insulators  do  not  affect  the  measurements  of  ioni¬ 
zation  currents.  The  leakage  current  betv/een  ti./o  electrodes  of  differing  poten¬ 
tial  can  be  elim.inated  by  using  a  split -insulator ,  one  section  of  v/hich  is  main¬ 
tained  at  the  potential  of  one  of  the  electrodes  .  Evidently,  no  leakage  path 
can  then  exist  across  this  equipotential  surface . 

The  three-electrode  ceramic-metal  chamber,  previously  employed  for  studies 
of  ion  source  rate  in  noble  gases,  was  modified  to  utilize  this  split-insulator 
principle  for  cesium  studies .  The  portion  of  the  center  electrode  that  was 
interior  to  the  tube  was  cut  away  to  leave  an  annular  Kovar  disc  which  was  level 
with  the  tvra  adjacent  ceramics  (Fig.  l)  . 
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ig  •  1  •  Ap-Cs  ion  ■tu])G  ochcinri.'ti''* 

If  a  potential  oource  is  applied  aero.-o  ai,  and  MB  Ir  ohorted  nlth  a  lead 
of  neglieible  resistance,  no  leakage  current  can  flov,  betueen  U  and  B.  Ibe  die 
tribution  of  electric  field  In  the  diode  between  U  and  B  will  not  be  signifi¬ 
cantly  altered  by  the  presence  of  M  because  the  ratio  of  the  gap  separation  to 
the  electrode  diameter  is  small  (  0.2). 

However,  in  order  to  measure  the  Ionization  current  in  the  diode  a  current 

meter  (A^^)  must  be  placed  in  series  with  U  and  B  and  adjacent  to  B  as  sho™  in 
Fig.  2. 


Fig.  2.  Circuit  that  v/ould  measure  ion  r-urrent  but  ar-  well  vrould 
measure  current  due  to  leakage  across  the  ceramic. 

Since  the  passage  of  ionization  current  through  produces  a  potential 

''mb'  ">  B  (0.1  V  for  full-scale  deflection  on  a  Ilexom  meter)  u 

feedback,  loop  mu.st  he  provided  in  order  to  maintain  thir-  potential  drop  as 
close  to  zero  as  possible.  Since  the  Hezem  is  already  provided  nth  a  potential 


I 
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compensator  in  the  form  of  a  -1  volt  full-rcalc  output ^  this  v/as  used  in  con¬ 
junction  v/ith  a  high-impedance  divider  and  a  Kintel  amplifier  (Model  ll4C)  to 
continuously  maintain  close  to  zero.  The  arrangement  is  shovm  in  Fig.  3. 
(The  maximuin  current  through  the  output  terminals  of  the  Hexem  meter  cannot 
exceed  120  [.xA  so  that  the  Hexem  alone  cannot  be  used  as  a  voltage  compensator 
for  ionisation  currents  >milliamps.) 


Fig.  3»  Circuit  used  to  nullify  influence  of  cesium  leakage 
across  ceramic  . 

The  current-voltage  characteristic  for  this  cesium  diode  is  then  given 
‘  "^UB  "'^MB  "  ^  ■ 

For  the  accurate  determination  of  ionization  currents  stringent  con¬ 

ditions  of  voltage  balance  are  often  required  betvrcen  M  and  B.  Experimental 
tests  with  a  cesium-filled  ceramic -metal  diode  indicated  that  the  leakage  resist¬ 
ance  between  M  and  B  can  be  as  low  as  -o  500  fl  .  (Most  often^  was 

found  to  be  — *  1000-2000  fL  for  liquid  cesium  temperatures  300°C.)  Now  if  the 
largest  error  that  can  be  tolerated  in  ionization  cui'rent  I  due  to  a  potential 
drop  is  1%,  then 
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®''mb  -  Rmb  I 

(for  I^g=500.a^I=10'*^A) 

=  5  X  10  ^  volts . 

Under  these  conditions  it  becomes  essential  to  use  separate  non-current 
carrying  leads  to  measure  Evidently  for  larger  values  of  the  con¬ 

ditions  of  voltage  balance  become  correspondingly  less  stringent . 

Ill .  RESULTS  OF  ARGON-CESIUM  TUBES  1  AliD  2 
A.  Tube  Design 

The  components  of  the  noble  gas  ion  tube(Pig::0,Sec.i^vrere  fjund  to ’,0  compatible 
with  cesium  up  to  bath  temperatures  around  250°C.  Operation  at  higher  tempera¬ 
tures  -  400“C  -  presented  problems  which  were  a  function  of  materials,  e  .g . 
choice  of  alumina,  and  particularly  of  pinch-off  tabulation.  The  latter  prob¬ 
lem  was  solved  by  employing  a  copper  rather  than  a  Kovar  pinch-off  as  previously 
used . 

mie  basic  tube  structure  was  that  of  KglO^Sec.A  with  the  center  electrode  cut 

out  to  allow  nulling  the  leakage  current  across  the  alumina  insulator  rings  due 

to  cesium  coverage  .  The  interelectrode  spacing  for  this  tube  was  6 .25  mm  or 

about  tv/ice  that  of  the  earlier  noble  gas  ion  tubes  because  of  the  removal  of 

the  center  electrode  and  the  use  of  the  original  ceramic  rings  as  insula¬ 
tors  . 

The  tube  was  mounted  in  an  oil  filled  four  foot  support  section  and  inside 
an  electrical  heater  assembly  as  seen  in  Fig.  4.  Electrical  heat  applied 
independently  to  three  heaters,  top  (above  the  tube),  middle  and  stem  (tubula- 
tion).  Chromel  alumel  thermocouples  were  mounted  on  the  uranium,  middle  and  bot¬ 
tom  electrodes  and  two  were  on  the  end  of  the  tabulation  (stem  thermocouples). 

B.  Argon-Cesium  Ion  Tube  1 

This  tube  was  processed  in  a  manner  similar  to  noble  gas  tubes  up  to  the 
backfilling  stage .  It  was  backfilled  with  argon  at  240  torr  and  sealed  from 
the  system  vath  the  cesium  chromate  pellets  attached.  After  seal-off  the  pel¬ 
lets  were  fired  and  the  cesium  was  transferred  into  the  tube  by  applying  heat 
to  the  pellet  vial.  The  copper  tabulation  was  then  pinched  off.  Due  to  the 
presence  of  the  argon,  the  cesium  did  not  transport  too  readily  and  considerable 
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Fig.  4.  Oven  for  argon-cesium  ion  tube,  v/ith 
tube  configuration  of  tubes  1  and  2 . 
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darken_^ng  of  "the  glasc  vial  wa;;  oboorved  . 
across  IhG  coramics  duo  bo  cGsiutn  coverage 


Tliorc  v/ar  no  por--eptihle  leakage 
.  It  v;ar  concluded  that  due  porrihlv 


to  difficulty  of  traneport  in  the  argon,  the  ceciura  condenred  on  the  tut,ulation 
before  entering  the  major  tube  volume.  Preouraably  cerium  could  be  on  the  tubu- 
lation  attached  to  the  sealed-off  ceramic  tube  ro  that  no  further  torts  were 
felt  necessary  to  'Ush  the  presence  of  cesium.  A  broakdovm  test  (as  a 
function  of  bath  temperature)  was  considered  but  not  run  duo  to  possible  damage 
to  the  uranium  electrodes.  Breakdovm  of  a  moekup  tube  did  change  as  a  function 
of  bath  temperature  Indicating  that  this  is  a  ccn-itlve  method  for  detecting 
the  presence  of  cesium  in  noble  gaser  . 

Typical  data  obtained  at  2  MV/  reactor  pov/er  is  shovm  in  Fig  .  5  . 


The  I-V  characteristic  appears  similar  to  pure  argon  data  taken  earlier, 
both  in  shape  and  in  magnitude.  Bie  asymmetry  in  the  I-V  data  is  consistent 
vath  the  fact  that  more  ions  are  generated  near  the  ur'ininm  electrode  and  sub¬ 
sequently  lead  to  a  higher  current  when  that  elor-trodo  is  negative.  The  cesium 
stem  temperature  for  these  data  was  2!30°C  .  However,  a  :  tem  temperature  varia¬ 
tion  from  135 °C  to  250°C  did  not  influence  appreciably  the  magnitude  of  the  cur¬ 
rent.  These  observations  suggest  that  the  tube  characteristic-  is  that  of  pure 
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argon  and  that  the  slight  influence  of  cesium  bath  temperature  on  tube  current 
may  be  due  to  cesium  at  much  less  than  its  equilibrium  vapor  pressure  . 

C  .  Argon -Cesium  Ion  Tube  2 

To  insure  the  presence  of  cesium  in  the  next  tube  it ‘was  processed  by  first 
firing  the  cesium  pellets  and  transporting  cesium  into  the  ion  tube  under  vacuum 
conditions.  Then  the  argon  w  .s  backfilled  at  240  torr  before  pinch  off.  Pres¬ 
ence  of  cesium  was  confirmed  by  a  dead  short  betv/een  electrodes  occurring  when 
the  cesium  was  shaken  into  the  ion  tube  from  the  tubulation.  This  tube  had  a 
3  mm  electrode  gap  by  the  use  of  thinner  ceramic  insulator  rings. 

In  order  to  operate  the  leakage  bucking  circuit,  a  ceramic  impedance  of 
several  hundred  ohms  is  needed.  Such  an  impedance  vras  not  obtained  by  thermally 
heating  the  tube  with  the  oven.  In  addition,  several  amps  at  low  voltage  v/ere 
passed  through  the  cesium  coating  in  order  to  burn  it  off.  Laboratory  tests  of 
this  indicated  an  impedance  of  hundreds  of  ohms  could  be  maintained  this  way 
and  the  tube  was  then  put  in  the  reactor . 

The  operation  of  this  tube  was  significantly  different  from  Tube  1.  I-V 
data  could  only  be  obtained  at  a  position  19  in.  out  fi’om  the  fully  inserted 
position  of  the  support  tube.  At  this  position  the  neutron  flux  is  approxi¬ 
mately  that  obtained  for  100  kw  reactor  power  v/ith  the  tube  fully  inserted  . 
Immediately  after  inserting  the  tube  fully,  the  leakage  resistance  of  the  ion 
tube  dropped  from  60  ohms  to  1  ohm  or  less.  Moving  the  tube  19  in.  out 
allowed  the  60  ohms  impedance  to  be  restored  by  'burning'  the  cesium  film  off 
the  ceramic  . 

The  current  magnitude  for  Tube  2  was  in  the  same  range  as  that  for  Tube  1 
for  equivalent  neutron  flux.  The  characteristic  followed  a  dependence 

out  to  20  volts  but  beyond  that  the  current  exhibited  an  increase  v/hich  may  be 
indicative  of  electronic  ionization.  Cesium  bath  temperature  did  not  influence 
the  current  magnitude  . 

These  observations  suggest  that  excess  cesium  v/as  present  in  the  tube  v/hich 
could  not  be  condensed  at  the  stem  cold  spot  due  to  the  poor  transport  proper¬ 
ties  of  cesium  in  argon.  Possibly  the  follovang  may  be  occurring.  At  burn-off 
from  the  ceramic,  the  cesium  condenses  onto  the  electrodes  of  the  diode.  Upon 
entering  a  high  radiation  field,  the  cesium  is  boiled  off  the  electrodes  by 
fission  fragments  and  recondensed  onto  the  ceramic ,  With  a  sufficiently 
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low  radiation  field,  the  evaporation  rate  of  cesium  from  the  ceramic  can  be 
maintained  greater  than  the  condensation  rate  .  T].ierefore  it  appears  that  data 
for  argon  with  cesium  at  some  equilibriiun  vapor  pressure  were  not  obtained . 

The  data  from  argon-cesium  ion  tubes  1  and  2  indicated  that  it  is  difficult 
to  maintain  cesium  at  an  equilibrium  vapor  pressure  in  240  torr  argon  in  a  high 
radiation  field  The  mean  free  path  of  neutral  cesium  in  argon  under  these  con¬ 
ditions  is  2x10“  cm  so  that  diffusion  times  will  be  long.  Transport  studies 
of  cesium  transport  in  argon  \ias  thus  undertaken . 


IV.  DIFFUSION  RATE  OF  CESIUM  VAPOR  IN  ARGON 


A .  Apparatus  and  Procedure 


The  apparatus  used  for  the  study  of  the  diffusion  of  cesium  in  argon  is 
shov/n  schematically  in  Fig.  6.  Tvro  molybdenum  electrodes  vrere  contained  in  a 
glass  chamber  filled  vdth  argon  v/hich  was  connected  to  a  liquid  cesium  reser¬ 
voir  by  a  glass  tube .  The  temperature  of  the  chamber  and  the  reservoir  section 
could  be  independently  controlled  by  means  of  separate  furnaces  (indicated  in 
Fig.  6  by  dotted  lines).  On  heating  the  liquid  cesium  to  a  constant  tempera¬ 
ture  T^  in  the  range  473  to  573 °K  the  diffusion  and  consequent  build-up  of 
cesium  vapor  in  the  chamber  (t^  ^  573 °K)  continuously  changed  the  value  of  the 
sparking  potential,  V^,  and  the  leakage  resistance,  R^,  between  the  electrodes. 
(R^  resulted  from  deposition  of  cesiujn  on  the  glass  between  the  electrodes.) 
When  V^  and  R^  became  constant  (at  the  same  time)  the  pressure  of  cesium  in  the 
chamber  was  considered  to  be  constant  and  at  the  value  of  the  vapor  pressure  of 
the  liquid  cesium  at  that  particular  temperature  . 


Tests  on  the  transport  of  cesium  in  vacuum,  and  in  argon  at  different 
pressures  (25  <  P  <  2J+0  torr),  indicated  that  for  constant  T  ,  R^  was  dependent 
only  on  p^  the  pressure  of  cesium  in  the  chamber;  i  .e .  P^  vras  not  a  function  of 
the  argon  pressure  (as  was  V^)  .  vath  this  knowledge  of  R^^  as  a  function  of  p  , 
the  follovang  procedure  could  be  adopted  for  transporting  cesium  from  the  reser¬ 
voir  to  the  chamber.  If  the  required  cesium  pressure  for  an  Ar-Cs  mixture  vritli 
Cs/Ar  of  10  was  vath  a  corresponding  bath  temperature  ^  (^  “5  523  °K  for 
p(Ar)  <  240  torr),  then  the  liquid  cesium  vras  heated  to  -o  573 °K  and  maintained 
at  this  value  until  R^  reached  a  value  corresponding  to  .  The  liquid  cesium 
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Fig.  6. 


oL'hemat  i  c 


^irrangemont  for  rtudier  of  diffiirion  of  oeriuiri  ’/ ipor  in  ai*jv.n 


temperature  was  then  quickly  lowered  to  the  value  O  ,  In  this  manner  the  dif- 

c 

fusion  rate  of  cesium  vfas  significantly  increased  over  that  which  vrould  have 
been  obtained  by  maintaining  the  lower  bath  temperature  of  O  . 

B.  Theory 

An  accurate  calculation  of  the  diffusion  coefficient  D  for  the  diffusion 
of  one  gas  (or  vapor)  into  another  gas  is  a  complicated  problem^^^  which  must 
take  into  account  the  nature  of  the  elastic  collisions  between  dissimilar  atoms 
and  also  atoms  of  the  same  kind  .  However,  a  simplified  expression  for  D  can  be 
obtained  if  the  concentration  n^  of  one  gas  is  alwaj/s  much  greater  than  the  con¬ 
centration  n^  of  the  other  gas,  and  if  the  momentum  loss  cross  section  Q^(a,c) 
for  collision  of  an  atom  of  each  kind  is  independent  of  relative  velocity.  Under 
these  conditions. 
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low  radiation  field,  the  evaporation  rate  of  cesium  from  the  ceramic  can  be 
maintained  greater  than  the  condensation  rate .  Tlierefore  it  appears  that  data 
for  argon  with  cesium  at  some  equilibrium  vapor  pressure  were  not  obtained  . 

The  data  from  argon-cesium  ion  tubes  1  and  2  indicated  that  it  is  difficult 
to  maintain  cesium  at  an  equilibrium  vapor  pressure  in  240  torn  argon  in  a  high 
radiation  f ield  The  mean  free  path  of  neutral  cesium  in  argon  under  these  con¬ 
ditions  is  2x10“  cm  so  that  diffusion  times  will  be  long.  Transport  studies 
of  cesium  transport  in  argon  was  thus  undertaken. 

IV.  DIFFUSION  RATE  OF  CESIUM  VAPOR  IN  ARGON 
A .  Apparatus  and  Procedure 

The  apparatus  used  for  the  study  of  the  diffusion  of  cesium  in  argon  is 
shown  schematically  in  Pig.  6.  I\to  molybdenum  electrodes  were  contained  in  a 
glass  chamber  filled  vdth  argon  which  ms  connected  to  a  liquid  cesium  reser¬ 
voir  by  a  glass  tube  .  The  temperature  of  the  chamber  and  the  reservoir  section 
could  be  independently  controlled  by  means  of  separate  furnaces  (indicated  in 
Fig.  6  by  dotted  lines).  On  heating  the  liquid  cesium  to  a  constant  tempera¬ 
ture  T^  in  the  range  473  to  573 °K  the  diffusion  and  consequent  build-up  of 
cesium  vapor  in  the  chamber  (T^  573 °K)  continuously  changed  the  value  of  the 

sparking  potential,  V^,  and  the  leakage  resistance,  R^,  between  the  electrodes. 
(Rj^  resulted  from  deposition  of  cesium  on  the  glass  between  the  electrodes.) 

When  V^  and  R^^  became  constant  (at  the  same  time)  the  pressure  of  cesiimi  in  the 
chamber  was  considered  to  be  constant  and  at  the  value  of  the  vapor  pressure  of 
the  liquid  cesium  at  that  particular  temperature  . 

Tests  on  the  transport  of  cesium  in  vacuum,  and  in  argon  at  different 
pressures  (25  p  <  240  torr),  indicated  that  for  constant  T  ,  R^  was  dependent 
only  on  p^  the  pressure  of  cesium  in  the  chamber;  i  .e  .  R^  was  not  a  function  of 
the  argon  pressure  (as  was  V  ).  With  this  knowledge  of  R  as  a  f motion  of  p  , 
the  follovang  procedure  could  be  adopted  for  transporting  cesium  from  the  reser¬ 
voir  to  the^chamber.  If  the  required  cesium  pressure  for  an  Ar-Cs  mixture  vrith 
Cs/Ar  of  10  was  with  a  corresponding  bath  temperature  ^  (^  523  °K  for 

p(Ar)^  240  torr),  then  the  liquid  cesium  was  heated  to  573°K  and  maintained 

at  this  value  until  R^  reached  a  value  corresponding  to  .  The  liquid  cesium 
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Fig.  6.  8ci.  matic  arrangement  for  studies  of  diffusion  of  cesium  vapor  in  argon 


temperature  was  then  quickly  lowered  to  the  value  O  ,  In  this  manner  the  dif- 

c 

fusion  rate  of  cesium  was  significantly  increased  over  that  which  would  have 

been  obtained  by  maintaining  the  lower  bath  temperature  of  0  . 

c 

B.  Theory 

An  accurate  calculation  of  the  diffusion  coefficient  D  for  the  diffusion 
of  one  gas  (or  vapor)  into  another  gas  is  a  complicated  problem^^^  which  must 
take  into  account  the  nature  of  the  elastic  collisions  between  dissimilar  atoms 
and  also  atoms  of  the  same  kind.  However,  a  simplified  expression  for  D  can  be 
obtained  if  the  concentration  n^  of  one  gas  is  always  much  greater  than  the  con¬ 
centration  n^  of  the  other  gas,  and  if  the  momentum  loss  cross  section  Q^(a,c) 
for  collision  of  an  atom  of  each  kind  is  independent  of  relative  velocity.  Under 
these  conditions. 
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D 


1/2 

jt  k  T(m  +  m  )  1 
a  c 


2mm 
a  c 


1 


u 


where  u  = 


8  k  T(m  +  m  ) 
a  c 


n  m  m 
a  c 


1/2 


(1) 


The  suffixes  a,c^  now  refer  to  argon  and  cesium,  respectively,  and  u  is  the 
average  relative  velocity  for  collisions  of  argon  and  cesium  atoms  . 

In  the  steady  state  the  current  density  of  cesium  atoms  due  to  diffu¬ 
sion  is 


Jc  =  -  . 


(2) 


It  is  assumed  that  the  time  taken  to  obtain  a  uniform  flow  rate  of  cesium 
through  the  0.45  cm  diameter  tube  (Fig.  6)  is  small  compared  vdth  the  time 
required  to  obtain  sufficient  cesium  in  the  chamber .  Since  the  temperature  of 
the  chamber  is  similar  to  the  temperature  of  the  cesium  bath  during  the  time 
for  diffusion  (see  above  procedure)  u  is  independent  of  position.  Tlierefore 
for  diffusion  along  tne  axis  of  the  tube, 

J  =  _  1  u  ^ 

c  3  dz  (3) 


where  is  the  current  of  cesium  atoms  per  second  passing  through  the  cros: 

sectional  area  cr  of  the  tube.  For  ^^c  =  0,  the  concentration  gradient 

dt  57“ 

is  independent  of  z  and 


n  (z)  =  n  (0)  +  z  ^^c 

dV  ■ 

Since  n^{L)  always  (L  is  the  length  of  the  tube), 


n 

c 


dn 

c 

dz 

at  z  =  0  is  given  by  the  value  of  n 


^  at  T  =  5T3°K. 
c  c 


(4) 
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Let  F  be  the  fraction  of  the  inride  curface  area  A  of  the  chamber  that  ir 
covered  vrith  cesium  v/hen  the  cesium  presrure  in  the  chamber  har  reached  an  equi¬ 
librium  value.  F(T  )  is  independent  of  p^  .  If  V  ir  the  volume  of  the  chamber, 

c  a 

then  the  total  number  of  cesium  atoms  that  have  t'  be  transported  from  the 
reservoir  is 


g 


V  + 


4  .F  A 

It  d^ 


v/here  d  is  the  atomic 

diameter  of 

ohe  cesium  atom . 

Under  these  simplified  conditions,  the  time  t  to  reach 

i  T  ) 
1  c 

n  (T  )V 
c  c 

+  4.F(T  )  .A 
c 

1  T 
'  g  ‘ 

1 

a 

‘8kT  (m  +m  )' 
g  c  a 

^/2  n^(573) 

0 

j 

\  % 

n  m  m 
c  a 

L 

equilibrium 


is 


Diffusion  times  in  the  chamber  are  considered  to  be  only  small  corrections  to  t 


C .  Results 


The  following  constants  are  used  for  evaluating  t  as  a  function  of  and 

n  for  the  particular  experimental  arrangements  shovm  in  Fig .  6 . 

^  -l6  /  o.,\  -1 


k 


m 


m 


g 


7t  m  m 
c  a 


V 

A 


(Massey  and  Burhop  ) 


n^(573) 


1.38  X  lO”"^^  erg  (‘’K) 

6.68  X  lO"""'^  g 
-22 

2.22  X  10  ^  g 
573  “K 


^4-1 

6.25  X  10  cm  sec 


3 


44  cm 


102  cm 

2 

0 .16  cm 

-14  2 

5.72  X  10  cm 

Q 

5.3*+  X  10~  cm 

2.8  X  10^^  atcans  cra“^ 

7.5  cm 
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Sub sti "tut/ ion  in  Ec[ .  (6)  yiolds 

t(sec)  ^  \[.03^  N^(T^)  +  2.09  F  (T^)] 


(6a) 


v/here  N  ,  N 


a- 


represent  the  number  of  argon,  oesium  atoms  cm"-'  in  units  of  10^*^ 


No  data  have  been  found  for  the  fractional  coverage  of  glass  by  cesium  as  a 
function  of  temperature  but  data  are  available  for  the  fractional  coverage  of 
tungsten  as  a  function  of  temperature^  ’ ^  for  T  -^420  to  2000°K.  These  data 
were  used  to  provide  a  measure  of  F(T^)  for  substitution  in  Eq.(r:a).  However, 
for  much  of  the  present  study,  the  fi,-t  term  in  Eq.(6a)  is  significantly 
larger  than  the  second  term  (provided  F(T_)  for  /^lass  is  of  the  same  order  of 
magnitude  as  F(T^)  for  tungsten)  so  thau  the  value  of  F(T^)  in  this  approximate 
calculation  is  unlikely  to  critically  affect  the  value  of  t .  Computations  of 
t  were  made  for  the  particular  conditions  under  which  the  experimental  values 
of  t  were  measured,  and  both  the  experimental  and  calculated  values  of  t,  in 
hours,  are  shovm  in  Table  1. 


TABLE  I 


O 

o 

N  T 
c  c 

(atoms  cm~^ 

X  10^^) 

F(T^) 

(fractional 

coverage) 

N 

(atoms  cm 

X  I0I6) 

Pa 

(torr  at  573 °K) 

t(exp) 
(Hours ) 

t(calc 

(Hours 

546 

1.50 

0.73 

42.2 

25 

0.9 

0.35 

550 

1 .65 

0.74 

42.2 

25 

1.2 

0.39 

544 

1.40 

0.73 

84.3 

50 

1.5 

0.66 

559 

2.12 

0.75 

169 

100 

2.5 

2.0 

48l 

0.21 

0.67 

169 

100 

0.5 

0.22 

503 

0.43 

0.69 

4  06 

240 

1.9 

1.0 

516 

0.65 

0.70 

4o6 

240 

2.0 

1.5 

The  agreement 

between  the  experimental  and 

calculated  value 

s  of  t  is 

con- 

siaered  to  be  satisfactory  in  view  of  the  approximations  made  in  the  calculation 
and  also  because  of  the  uncertainty  in  the  experimental  values  of  t  due  to  the 
difficulty  of  accurately  and  quickly  changing  temperature  .  In  most  of  the  exam¬ 
ples  shown  in  Table  I  the  ratio  of  cesium-argon  N^)  is  less  than  lO"^ . 

If  this  ratio  is  increased  to  the  desired  value  by  decreasing  N  .T  ,  then  the 
value  of  t  is  correspondingly  decreased. 
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It  appears^  therefore,  that  the  present  calculations  can  yield  useful  values 
of  the  times  required  to  obtain  an  equilibrium  cesium  pressure  in  an  ionization 
chamber  due  to  diffusion  of  cesium  vapor  from  a  cesium  reservoir,  and  that  such 
times  should  be  ^  2  hours  (’vith  the  apparatus  considered  here)  for  a  mixture  of 
argon-cesium  vsith  Cs^  Ar  of  10  ^  for  argon  pressures  ^  240  torr  . 

V.  RESULTS  OF  ARGON-CESIUM  TUBE  3 
A .  Tube  Design 

The  cesium  transport  studies  indicated  that  a  relatively  large  diameter 
tabulation  was  required  betv/een  the  cesium  reservoir  and  the  ion  chamber  shovm 
in  Fig .  4 .  In  order  to  ensure  that  the  increased  opening  to  the  ion  chamber 
(1/^  •  o ‘d  .  copper  tubing)  did  not  interfere  vrlth  the  measurements  of  the 

diode  characteristic,  a  Kovar  electrode  was  supported  in  the  manner  shovm  in 
Fig .  7  so  that  ion  collection  on  a  planar  Kovar  electrode  could  still  be  made  . 

The  spacing  for  this  ion  tube  v/as  set  at  5  mm  •  By  separating  the  ceramic  rings 
with  the  Kovar  ring  r-hovm  in  Fig .  7  sufficient  area  v;as  provided  so  that  cesium 
could  be  transported  readily  from  the  cesium  reservoir  to  the  volume  between 
the  uranium  and  Kovar  electrodes  . 

The  circuitry  and  heater  were  similar  to  those  discussed  in  Sections  II 
and  III .  However,  for  this  tube  considerable  data  were  obtained  v/ithout  the 
need  to  use  the  nulling  circuit  since  the  leakage  currents  across  the  ceramics 
were  much  less  than  the  measured  ion  currents . 

B .  Experimental  Data 

This  ion  tube  operated  satisfactorily  in  the  reactor  and  the  data  obtained 
were  reproducible  once  the  cesium  vapor  equilibrium  condition  vras  reached  .  In 
order  to  determine  the  influence  of  the  rather  large  area  of  the  middle  elec¬ 
trode,  data  were  taken  vrLth  the  middle  electrode  connected  to  the  uranium  elec¬ 
trode  and  then  to  the  Kovar  electrode .  These  data  indicated  that  v/hen  ion 
current  is  collected  on  a  single  electrode,  either  uranium  or  Kovar,  the  cur¬ 
rent  is  well  behaved.  That  is,  the  current  magnitude  was  as  might  be  expected 
and  under  most  conditions  the  current  follov/ed  a  relationship.  There  was  a 

significant  displacement  of  the  zero  current  from  the  origin  of  the  voltage 
scale  however  (0.5  volts  with  Kovar  negative  and  uranium  positive)  which  is  not 


14 


completely  understood.  (The  noble  ga::  data  alv/ays  went  througli  zero  vollr. ) 
This  may  be  associated  v;ith  some  small  thennionii''  current  which  appeaj’ed  to 
influence  the  ion  current  only  v/hen  the  tube  was  fully  inserted  into  the  core 
at  2  MV/  po’wer  level .  In  thi  s  case  the  i-V  characteristic  v/as  very  as^/mmetric — 
the  uranium  negative  side  of  the  characteristic  shovring  100  timer;  more  cuirrent 
than  the  other.  Plotted  on  a  logarithmic  scale  the  characteristic  is.  s.hovm  in 
Fig  .  8 . 


In  order  to  make  an  estimate  of  surface  conditions  that  could  give  rise 
to  such  electron  currents,  the  follov/ing  assumptions  were  made. 

1.  The  slope  of  the  retarding  field  curve  could  be  used  to  obtain 
the  uraniiam  temperature ; 

2.  The  first  current  saturation  at  25  ma  corresponds  to  saturated 
emission,  and 

3.  The  knee  of  the  curve  gives  the  difference  in  v/ork  function  between 
the  uranium  and  Kovar  surfaces . 

From  the  slope,  T^=850°K  v/hich  for  a  current  of  25  ma  corresponds  to  a  surface 
of  work  function  1 .7  ev .  This  implies  that  the  cesium  coating  on  the  uranium 
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provide  oUch  a  lov/  work  function  and  aloo  that  a  larirre  temperature  gradient 
exists  between  the  point  where  the  temperature  was  measured  (the  edge  of  the 
disc)  to  the  uranium  surface  inside  the  tube.  Such  a  possibility  is  not  mirea- 
sonable.  Now  the  collector  work  function  (Kovar)  can  be  determined  from  tie 
absolute  magnitude  of  the  retarded  current  at  a  particular  voltage;  this  gave 
a  value  of  1 .8  ev  -  a  number  which  appears  to  be  reasonable  for  cesium  on  Kovar. 
The  knee  of  the  curve  is  roughly  at  0.1  volts  which  corresponds  well  to  the  dif¬ 
ference  of  the  computed  work  functions  for  uranium  and  Kcvar.  Thus  it  -appears 
that  cesium  coverage  of  the  electrodes,  along  vath  a  temperature  higher  than 
that  expected  on  the  uranium  surface,  can  explain  the  current  asymmetry  attribu¬ 
ted  to  thermionic  emission  from  uranium. 

In  order  to  avoid  this  electron  emission  effect,  the  tube  v/as  operated  at 
10  inches  above  the  core  v/here  the  uranium  v/as  at  a  much  lower  temperature. 

The  analysis  of  this  v/ell  behaved  data  follov/s . 


C.  Pit  of  Theory  to  Experimental  Data 

Unlike  the  case  for  the  neon-argon  mixture  a  theoretical  current -voltage 
curve  could  not  be  computed  for  the  argon-cesium  mixtures  since  W,  the  energy 
to  create  an  ion  pair,  has  not  been  reported  in  the  literature  for  this  gas 
mixture .  Instead  a  typical  inpile  experimental  run  was  selected  and  a  one 

parameter  fit  was  made  to  the  data  by  adjusting  W  in  the  code  solution.  The 
fit  is  shovm  in  Fig .  9  . 

The  inpile  run  selected  for  comparison  v/as  made  v/ith  the  reactor  at  2  MV/ 
power  but  the  ionization  tube  v;as  withdraw  10  inches  from  the  midplane  of  the 
core.  The  cesium  bath  temperature  was  150"C  which  for  an  argon  pressure  of 
240  torr  gave  a  cesium  to  argon  ratio  of  Cs/Ar=2 .SxlO'^ .  For  this  position  the 
thermal  neutron  flux  had  been  measured  as  6.5x10^^  cm"^  sec"^  vath  gold  foil 
counting  techniques.  The  experimental  curve  (Pig .  9)  is  similar  in  shape  to 
the  neon-argon  I-V  curves;  that  is  the  current  follov/s  a  dependence  on 

voltage  out  to  about  +  25  volts  and  then  the  current  increases  rapidly  v/ith 
voltage  again  presumably  due  to  electron  impact  ionization  in  the  sheath .  The 
code  curve  v/as  adjusted  to  fit  the  experimental  data  in  the  negative  quadrant 
(-25  volts  ^  V^O  volts),  v/here  the  positive  ion  sheath  is  adjacent  to  the 
Kovar  electrode.  The  fission  fragment  ranges  for  argon  v/ere  used  and  the  diode 
spacing  v/as  O.504  cm.  As  in  the  case  for  the  noble  gas  code  runs  (Section  A)  the 
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Gas  Pressure  •  240  Torr 


Neutron  Flux  •  6.5  x  10*^  sec'*  cm'^ 


Fig.  9.  l-v  data  for  argon-cesium  tube  3  compared  to  theory  using 
as  an  adjustable  parameter. 

ion  collection  area  for  V  0  v/as  taken  as  the  total  Kovar  electrode  area(Trp2^) 
while  for  V  *-0  the  ion  collection  area  was  taken  as  just  the  uranium  area(np^^)  . 
The  code  curve  gave  the  fit  shovm  vath  a  value  of  ^^=1^+ .0  ev/ion  pair  or  ^  =14.84 
and  t72=15.54  ev/ion  pair.  Also  shovm  for  comparison  is  the  theoretical  curve  for 
pure  argon  (W^=26 .4  ev/ion  pair)  for  these  same  conditions. 

It  is  seen  that  the  current  for  Ar-Cs  is  greater  than  that  for  Ar  alone 
which  suggests  a  higher  ion  generation  rate  in  this  mixed  gas  system.  This  may 
arise  from  the  additional  ionization  of  cesium  from  the  argon  metastable  states  . 
On  the  other  hand,  Cs'*’  ions  may  diffuse  from  the  b  ilk  of  the  plasma  into  the 
sheath(as  discussed  in  Section  A  for  the  Ne-Ar  data)  and  thereby  gi-e  a  current 
higher  than  that  due  to  the  total  ions  generated  rate  within  the  sheath.  More 
detailed  investigations  are  under  v/ay  to  resolve  these  unknovms  . 
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SECTION  F 


SECONDARY  EIFCTRON  YIELD  FROM  FISSION  FRAGMENTS'^- 


ABSTRACT 

The  secondary  electron  yield  from  fission  fragments  has  been  measured 
directly  using  an  ultra  high  vacuum  diode  containing  uranium  v;hich  was  oper¬ 
ated  in  a  nuclear  reactor.  The  current  measured  in  this  diode  consisted  of 
contributions  from  the  positively  charged  fission  fragments,  secondary  elec¬ 
tron  emission  induced  by  fission  fragments  and  photoelectrons  generated  by 
gamma  radiation.  The  secondary  electron  yield, was  determined  from  the 
secondary  electron  current  and  the  computed  fission  fragment  emission  rate  . 

For  a  thick  (>one  range)  uranium-nickel  (5.7  v/t  %  nickel)  alloy  the  value 
of  A  is  207  +  10 . 


Secondary  Electron  Yield 

Measurements  have  been  reported on  the  secondary  electron  yield 
accompanying  the  emission  of  fission  fragments  from  uranium  dioxide  films  of 
thickness  0.1  to  3  microns.  The  average  range  for  the  light  and  heavy  frag¬ 
ments  in  uranium  dioxide  is  10.6  and  8.2  microns,  respectively.  The  elec¬ 
tron  yield  decreased  vrLth  increasing  film  thickness  from  about  57O  to  3OO 
electrons  per  fragment.  This  note  discusses  secondary  electron  yield  for  a 
film  of  a  uranium  metal  alloy  that  was  several  times  the  fission  fragment 
range  in  thickness  and  therefore  equivalent  to  a  yield  for  a  film  of  one 
range  thickness. 

The  secondary  yield  v/as  measured  with  an  ultra  high  vacuum  (lO”^  torr) 
ceramic-metal  tube  which  consisted  of  three  Kovar  electrodes,  each  of  2.5  cm 
diameter,  forming  two  diodes.  This  tube  v/as  developed  primarily  for  inpile 
studies  on  the  ionization  of  noble  gases  by  fission  fragments . The  uranium 
v^as  brazed  to  one  electrode  in  the  form  of  a  uranium-nickel  alloy  (5.7  wt  % 
nickel)  v/hich  was  I5  .5  microns  thick  and  1 .9  cm  in  diameter.  The  interelec¬ 
trode  spacing  in  the  diode  containing  uranium  was  O.30  cm  whereas  the  spacing 
in  the  diode  without  uranium  was  O.33  cm. 

•^Submitted  to  the  Journal  of  Applied  Physics  .  ^  - - - - 
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The  conventional  current  I  measured  in  the  uranium  diode  can  consist  of: 

(1)  Positively  charged  fission  fragments,  (ip; 

(2)  Secondary  electrons  emitted  from  the  uranium  by  fission  fragments 
leaving  the  uranium  electrode  (l^j 

(3)  Secondary  electrons  emitted  from  the  Kovar  by  fission  fragments 
arriving  at  the  Kovar  electrode  (l‘); 

(4)  Electrons  generated  by  gamma  collisions  (photo  and  Compton  effects) 
\d.thin  the  tube  walls,  (ipj  and 

(5)  Tertiary  electrons  generated  by  impact  of  electrons  in  processes 

2,  3^  ^  • 

The  diode  that  did  not  contain  uranium  (gamma  chamber)  was  used  to  study 
the  magnitude  of  processes  (Ij-)  and  (5)  at  a  gamma  flux  of  I.7  x  lo"^  Rad  hr"^. 
Figure  1  shows  the  current -voltage  characteristic  obtained  with  this  diode  where 
the  inserts  represent  the  particle  flows  at  large  negative  and  positive  voltages . 
The  current  is  small,  and  the  lack  of  saturation  at  high  voltage  indicates  a 
broad  electron  energy  distribution  characteristic  of  gamma -produced  electrons. 

The  photoelectrons  which  are  not  completely  retarded  by  a  high  retarding  voltage 
are  indicated  in  Fig .  1  by  the  dashed  curves .  Tbe  small  inflection  in  the  cur¬ 
rent  near  zero  voltage  suggests  the  presence  of  a  few  tertiary  electrons  (pro¬ 
cess  5)  and  the  asymmetry  in  the  curve  suggests  a  difference  in  photoemission 
from  the  two  surfaces.  However,  the  purpose  of  using  the  gamma  chamber  was  to 
estimate  the  gamma-induced  contribution  to  the  current  measured  with  the  uranium 
diode,  and  as  indicated  below,  this  gamma  contribution  is  not  significant. 

Figure  2  shows  the  results  for  the  diode  chamber  containing  uranium.  It 

may  be  seen  that  currents  measured  with  this  chamber  at  a  neutron  flux  of 
12  -1  -2 

5  *3  X  10  sec  cm  are  about  two  orders  of  magnitude  larger  than  the  currents 
measured  simultaneously  in  the  gamma  chamber  (Fig .  l) .  The  gamma  induced  cur¬ 
rent  will  therefore  be  a  small  correction  in  the  uranium  chamber.  At  large  posi¬ 
tive  potential  (Fig.  2),  we  see  that  =  I  +  The  secondary  electron 

yield  A  Per  fission  fragment  for  the  uranium  is  given  by  I^/eN  ,  where  N  is  the 
fission  fragment  particle  emission  rate.  The  current  at  large  negative  poten¬ 
tial  could  be  used  to  determine  secondary  electron  yield  per  fission  fragment 
for  Kovar;  however,  for  that  case  a  correction  would  have  to  be  made  for  the 
energy  and  angular  distribution  of  the  fission  fragments  that  bombard  that  elec¬ 
trode.  Note  that  tertiary  electrons  are  always  suppressed  at  large  potentials. 
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In  addition  the  saturation  characteristic  of  Fig .  2  indicates  that  the  secondary 

electrons  generated  by  fission  fragments  have  an  energy’’  distribution  lying  almost 

entirely  below  20  ev;  thus  tertiaries  will  be  few  in  any  event. 

+ 

The  current  from  fission  fragments  is  computed  from 

ij:  =  A<t  e(Z^R^  + 


where  A,  4),  represent  respectively  the  uranium  area,  the  neutron  flux,  and 
the  macroscopic  fission  cross  section;  eZ^  p  and  R  are  the  initial  charge 
and  range  of  each  fragment .  We  obtain  =  1 .7  x  10-'^  amp  for  the  data  of 
Fig,  2,  about  5/°  of  the  measured  saturation  current. 

The  fission  fragment  emission  rate  was  computed  from  N^.  =  A4>  +  ^2^^^ 

which  is  derived  in  the  Appendix  along  \d.th  Eq.  1.  The  value  of  current  obtained 
at  50  volts  is  used  to  calculate  the  secondary  electron  emission  coefficient . 
Table  1  lists  the  values  computed  for  two  different  flux  levels . 


TABLE  1.  Secondary  Electron  Yield  per  Fission  Fragment, A 

(Electrons/Fission  Fragment) 


Neutron  Flux 

(sec~^  cm.~^) 

5  .3  X  10^^ 

12 

5.3  X  10 


Uranium 
Electrode, A 
212  +  10 
202  -I-  10 


Anno^^^  has  indicated  an  extrapolated  value  for  A  about  100  for  a  one 
range  film  of  UO^ .  His  measurements  ext.ended  from  O.Y^o  to  28 7o  of  the  fission 
fragment  range  in  UO^ . 

Tliere  is  considerable  similarity  in  the  energy  distributions  of  electrons 


emitted  from  the  uranium  and  Kovar  surfaces,  because  the  current -voltage  char¬ 
acteristic  is  symmetric  about  the  current  at  sero  voltage .  It  can  be  readily 
checked  that  space  charge  effects  are  negligible  for  the  experimental  condi¬ 
tions  reported  here . 


APPENDIX 

The  fission  fragment  emission  rate  N  is  computed  assuming  the  fragments 
are  generated  uniformly  and  are  emitted  isotropically  at  the  point  of  genera¬ 
tion.  The  fragments  born  within  an  increment  dy  at  a  depth  y  below  the  uran¬ 
ium  surface,  and  moving  at  an  angle  0  with  respect  to  the  normal,  will 
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contribute  to  the 
ness) . 


emission  when  0^0 

max 


cos"^(y/R)  and  y 


R  (pne  range  thick- 


(1/2)  sin  ^  d  O  dy 


(2) 


o  0 

This  yields  N  =  A  <ti  R/4.  Since  two  groups  of  fragments  are  generated  per 
fission,  the  total  fission  fragment  emission  rate  is  N,  =  A  <i)  2  (r  +  R  )/4. 

mu  .  +  t  f'  1  2  ' 

The  fission  fragment  current  I  from  the  uranium  surface  is  obtained  by 
including  in  Eq.  (2)  an  expression  which  takes  into  account  the  linear  loss  of 
fragment  charge  with  distance. The  electronic  charge  on  a  fragment  escap¬ 
ing  from  the  uranium  surface  at  angle  O  from  depth  y  is  Z(l  -  y/R  cos  O)  where 
Z  is  the  initial  electronic  charge  of  the  fragment.  Thus, 


I*"  =  A  <t.  e  Z  r  r  (1  .  y/R  cos  0)(l/2)sin  ©dOdy,  (3) 

0  “^o 

which  yields  I  =  A  <ti  e  Z  R/8.  For  the  two  groups  of  fragments,  the  total 
fission  fragment  current  is  that  given  in  Eq..  (l)  of  the  terb . 
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